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ABSTRACT 

We use VLT/FLAMES intermediate resolution (R~6500) spectra of individual red gi- 
ant branch stars in the near- infrared Call triplet (CaT) region to investigate the wide- 
area metallicity properties and internal kinematics of the Sextans dwarf spheroidal 
galaxy (dSph). Our final sample consists of 174 probable members of Sextans with 
accurate line-of-sight velocities (± 2 km s _1 ) and CaT [Fe/H] measurements (± 0.2 
dex). We use the Mg I line at 8806.8 A as an empirical discriminator for distin- 
guishing between probable members of the dSph (giant stars) and probable Galactic 
contaminants (dwarf stars). 

Sextans shows a similar chemo-dynamical behaviour to other Milky Way dSphs, 
with its central regions being more metal rich than the outer parts and with the more 
metal-rich stars displaying colder kinematics than the more metal-poor stars. 

Hints of a velocity gradient are found along the projected major axis and along 
an axis at P.A.= 191°, however a larger and more spatially extended sample may be 
necessary to pin down the amplitude and direction of this gradient. 

We detect a cold kinematic substructure at the centre of Sextans, consistent with 
being the remnant of a disrupted very metal poor stellar cluster. 

We derive the most extended line-of-sight velocity dispersion profile for Sextans, 
out to a projected radius of 1.6 deg. From Jeans modelling of the observed line-of-sight 
velocity dispersion profile we find that this is consistent with both a cored dark matter 
halo with large core radius and cuspy halo with low concentration. The mass within 
the last measured point is in the range 2-4 x 10 s M©, giving very large mass-to-light 
ratios, from 460 to 920 (M/L)y )S . 

Key words: stars: abundances - galaxies: kinematics and dynamics - galaxies: dwarf 
- Local Group - Galaxies: individual: Sextans dSph - dark matter 



1 INTRODUCTION 

The Sextans dwarf spheroidal galaxy (dSph) is a satel- 
lite of the Milky Way ( heliocentric distan ce of 86 kpc 



lite or tne Milky w ay ( Heliocentric distance ol 80 Kpc 
fromllrwin et al.lll990l and lMateo et al.ll 19951 : 95.5 kpc from 
iLee et al.1 120031) . Sextans was discovere d relatively recentl y 
by automatic scanning of APM plates bv llrwin et ail (|l990t ). 



* Based on FLAMES observations collected at the ESO, proposal 
171.B-0588. 

f Corresponding author. E-mail: gbattagl@eso.org 



This is because the Sextans dSph (hereafter, Sextans) is one 
of the most diffuse and faint dSphs, with a central surface 
brightness Eo = 18.2 ± 0.5 mag arcmin -2 and luminos - 
ity L = 4.1 ± 1.9 x 10 5 L^ lllrwin fc Hatzidimitrioulll995l '). 
This, together wit h its locatio n on the sky (I = 243.5°, 
b = +42.3° from iMated 119981 ). results in a considerable 
amount of contamination from Milky Way stars. This makes 
it particularly difficult to explore the 2D structure of Sex- 
tans and search for possible a symmetries or other features 
ijlrwin fc Hatzidimitrioulll995h . 

A further difficulty in the study of Sextans is given by 
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its large extent on the sky (r core = 16.6±1.2 arcmin; r t idai = 
160.0±50 arcmin from I Irwin fc Hat zidimitrioulli995l ). which 
is not well suited to the field of view of most current wide 
area imagers and spectrographs. Several imaging studies of 
the central r egion (at most the i nner 40 arcmin x 40 arcmin) 
of Sextans jMateo et al.lll99ll. [l995l; iBellazzini et al.l [200ll : 
lHarbeck et al.ll200ll ; I Lee et all 120031 ) have shown that the 



from 504 members gives a mass of 5.4xl0 7 M0within the 
last measured point at ~65 arcmin (equivalent to ~4 core 
radii, ~0 .4 tidal radii), i .e. a m ass-to-light (M/L) ratio of 130 
(M/ LWlStrigari et all J2007h used a sample of 294 members 
from lWalker et alj 112006) and estimated a mass of 0.9_ n \ x 



majority of its stars are ancient (>10 Gyr old), as evidenced 
by the presence of horizontal branch (HB) and RRLyrae 
stars. A significant population of stars have also been found 
on the main sequence (MS), above the oldest turn-off. It 
was sugges ted that the se may be MS stars as young as 2 Gyr 
(e.g. lMateo et al.lll99l[ ). H owever, deep Suprime -Cam V and 
I photometry at Subaru l|Okamoto et alj|2008h shows that 
their spatial distribution is very similar to old (>10 Gyr) 
MS stars below the oldest turn-off, indicating that these are 
likely to be old blue stragglers (BS) stars (>1 Gyr) rather 
than young MS stars (see also I Lee et al.ll2003T ). 

In common with other dSphs, Sextans presents spatial 
variations of its stellar population mix, with the red horizon- 
tal branch (RHB) stars more centrally concentrated than the 
blue horizontal branch (BHB) star s IBellazzini et al.ll200ll ; 
lHarbeck et al ] |200ll ; [Lee et al. 2003; Kle vna et al.ll2004h . By 
estimating the luminosity difference between the HB and 
the main-se quence turn-off o ut to the nominal tidal radius 
of Sextans, lOkamoto et al.l (|2008j ) estimate an age differ- 
ence of at least 3 Gyr between the outer parts (with age 
~ 14 Gyr) and the central parts (with age ~ 10 Gyr). The 
younger stars are more centrally concentrated than the older 
ones, consistent with the different distribution of RHB and 
BHB stars. The other populations, red giant branch (RGB), 
BS and MS stars, have a spatial distribution intermediate 
to the one of RHB and BHB stars: this indicates that their 
spatial distribution can be thought of as a combination of 
the stellar distribution of RHB and BHB stars, i.e. RGB, 
BS and MS are likely to represent both components. 

No determinations have yet been made of the large scale 
metallicity properties in Sextans. Spectroscopic estimates 
of the mean metallicity ([Fe/H]) of Sextans are available 
from Call triplet (CaT) observations of individual red giant 
branch (RG B) stars in the centra l regions (<20 arcmin) of 
the galaxy JPa Costa et al.lll99ll . <[Fe/H]>= -1.7 ± 0.25 
dex, 6 stars; Isuntzeff et al.lll993l -2.05±0.04 dex, 43 stars), 
which also revealed an internal [Fe/H] spread. High resolu- 
tio n spectroscopic obser vations of 5 RGB stars in Sextans 
bv lShetrone et all (|200lf ) confirmed both the average [Fe/H] 
found in the CaT studies, and the large range of observed 
metallicities (A [Fe/H] = 1.4 dex). T he derivation of the 
mean [Fe/H] and spread in Sextans bv lGeisler fc Saraiedinil 
jl996t ) from Washington photometry is also in agr e ement 
with the spectroscopic studies. Recently lAoki et all (|2009T ) 
and Tafelmayer et al. (2010) follow-up at high spectral reso- 
lution 6 and 2 stars in Sextans, respectively, confirming the 
presence of stars with [Fe/H] down to —3.1. 

Studies on the internal kinematics of Sextans have 
shown that, in common with the other dSphs, this galaxy 
exhibits a much larger velocity dispersion than what is 
expected from the gravitational contribution of its lumi- 
nous matter, hinting at the presence of large amounts of 
dark matter jS untzeff et al. 199<§^_ Hargreaves et all 1 19941 ; 
iKlevna et al.ll2004l:| Walker et alj|2006h . The latest estimate 
(Walker et al. 2007) based on a velocity dispersion profile 



10 7 M within 0.6 kpc (13^ x 10 7 M Q extrapolated to 4 
kpc), which corresponds to an M/L= 260 (M/L)©. Although 
the mass estimates are not directly comparable as they refer 
to different distances from the Sextans centre, they all point 
to mass-to-light ratios so large as to require considerable 
amounts of dark matter, if the hypothesis of dynamical equi- 
librium holds. These spectroscopic samples do not reach the 
nominal tidal radius and therefore the masses determined 
are likely to be lower limits. 

With its unusually large extent on the sky and large 
tidal-to-core radius ratio (rtidai/fcorc ~ 9.6 as compared to 
ratios between 3 and 6 for most of the other classical dSphs, 
see Irwin & Hatzidimitriou 1995), Sextans seems a natural 
candidate for tidal disruption. However, no clear signs of dis- 
ruption such as tidal tails and S-shaped contours have been 
found in Sextans. This may however just be a consequence 
of the observational challenge of detecting such signs in a 
low surface brightness and heavily contaminated object like 
Sextans, or that such signs are not necessarily expected in 
tidally d isrupted galaxies, a s some N-body simulations show 
(see e.g. iMunoz et al.l [20081 ). A small velocity gradient has 
been detected ap proximately along th e projected minor axis 
of this galaxy by I Walker et al. (2008). However, given that 
no proper motions have been directly measured for Sextans, 
it is unclear whether this gradient is a sign of tid al disrup- 
tion a s predicted, for example, in the models o f lOh et al.l 
l|l995l ), due to geometrical effects |Strigarill20i(ih or an in- 
trinsic gradient. 

In this paper we use a spectroscopic sample of line-of- 
sight velocities and [Fe/H] measurements for 174 RGB stars 
probable members of Sextans, obtained at VLT/FLAMES, 
to analyse both the wide-field metallicity properties and the 
internal kinematics of this galaxy. Our sample extends out 
to 6 core radii (0.7 times the nominal tidal radius). The 
structure of the paper is the following. In Sect. [5] we de- 
scribe the observations, the data reduction technique and 
the adopted metallicity scale. In Sect. [3] we deal with the is- 
sue of separating Galactic contaminants from Sextans mem- 
bers and we use the Mg I line as an empirical discriminator 
between dwarf and giant stars. In Sect. [3] we present the 
results from our survey regarding the wide-area metallicity 
and kinematic properties of Sextans as well as the presence 
of kinematic substructures. Finally, we perform a mass de- 
termination in Sect. [S] and present discussion and conclu- 
sions in Sects. [6] and [7] 



2 OBSERVATIONS AND DATA REDUCTION 

The spectroscopic observations were carried out at 
VLT/FLAMES between 2003-2004, with the exception of 
one pointing which was re-observed in May 2008, during 
the commissioning of the upgraded GIRAFFE spectrograph 
CCD. 

For the selection of our spectroscopic targets we used 
photometric data from INT/WFC covering approximately 
an area of 2 degx 2 deg centred on the coordinates of the op- 
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tical centre galaxy and ESO/WFI photometry for 12 point- 
ings spread over the outer regions of the system. The imag- 
ing data were reduced following standard p rocedures, for 
details see references in lBattaglia et all (|2006l ). In this work 
these photometric data are used exclusively to select spec- 
troscopic targets. 

We selected targets classified as stellar in our photom- 
etry and with a position on the colour-magnitude diagram 
(CMD) consistent with an RGB star, but with a wide colour 
range to avoid biasing our sample in age or in metallicity. 

We used VLT/FLAMES feeding the GIRAFFE spectro- 
graph in Medusa mode, that allows the simultaneous alloca- 
tion of 132 fi bres (including sky f ibres) over a 25' diameter 
field of view i|Pasquini et alj [2002l. We used the GIRAFFE 
low resolution grating LR8 (resolving power R~ 6500) , cov- 
ering the wavelength range from 8206 A to 9400 A. This 
allows the measurement of equivalent widths (EW) from the 
near-infrared CaT lines at 8498, 8542, and 8662 A and also 
enables the derivation of velocities accurate to a few km/s. 
The spectra used in this work are from 16 different fields in 
Sextans obtained with this set-up (see Tab.[2]for the observa- 
tions journal). These data were reduced using the GIRAFFE 
pipeline (Geneva Observatory; Blecha et al 120031 ) . 

The sky-subtraction and extraction of the velocities and 
EWs of CaT lines were carried out using o ur own software 
develo ped by M.Irwin. For the details see Battaeli a et al.l 
(2008). Here we just note that we estimate the CaT EW in 
two ways. The first consists in simply summing the flux con- 
tained in a region 15A wide centred on each CaT line (here- 
after, integral fit). To derive the second estimate we fit in- 
dividual unconstrained Gaussian functions to each CaT line 
over the same wavelength region (hereafter, Gaussian fit). 
The combined EWs for CaT lines #2 and #3 (As542, -W62) 
for both the integral and Gaussian fits are then compared 
and used to compute an overall correction to the Gaussian 
fit. This is necessary since the observed CaT lines have non- 
Gaussian wings which are progressively more visible as the 
EW increases. 

Beside the EW of the CaT lines, we also measure the 
EW of a Mglline at 8806.8 A . The EW of this line is esti- 
mated using the values from the integrated fit, summing the 
flux contained in a region 6A wide centred on the Mg I line. 
The EW of this line will be used as a discriminant between 
Sextans stars and Milky Way contaminants as described in 
Sect. [M] 

We observed a total of 1294 targets, which include 103 
stars with double measurements, 70 with 3 measurements, 
5 with 4 measurements. The number of distinct targets is 
1036. 

As described in our previous studies using the same 
instrument and set-up we ensured a high reliability for the 
data used in our analysis by excluding those stars with S /N 
per A < 10 and vel ocity errors larger than 5 km s _1 (e.g. 
iBattaglia et al.l2008l ). In order to ensure the reliability of the 
EW measurements we also required that the summed EW of 
the two strongest CaT lines should not differ by more than 
2 A between the determination from the integral fit and the 
Gaussian fit. 

Figure [T] shows the distribution of differences in velocity 
for the stars that have double measurements (103 objects), 
both without any selection criteria and adopting the above 
selection criteria for each measurement (59 objects). Fig- 



ure [2] refers to the distribution of differences in summed EW 
from Gaussian fit of CaT lines and integrated flux for the 
Mg I line. The determination of the summed CaT EW from 
the integral and Gaussian estimators yield quite similar dis- 
tributions, and therefore, for consistency with our previous 
work, hereafter we use the summed CaT EW derived from 
the Gaussian estimator. For the EW of the Mglline the in- 
tegral estimator yields slightly less noisy measurements and 
therefore we will use this for the Mglline. The error in the 
summed CaT EW for the individual measurements is given 
by ffsw =6/(S/N), while the error in the MglEW is well 
represented by 2.8/(S/N). 

In the following, when calculating the velocity, summed 
CaT EW and all the other properties of the stars with re- 
peated measurements, we use only the individual measure- 
ments which meet our S/N and velocity errors criteria (be- 
fore combining them). 

The final sample was carefully checked to weed out any 
spurious objects (e.g. broken fibres, background galaxies, 
foreground stars, etc.). We removed 4 objects because they 
were observed with broken fibres; we found 3 background 
galaxies and removed 9 objects because the continuum shape 
or the presence of very broad absorption line was not con- 
sistent with what is expected for RGB stars. 

Our final sample of acceptable measurements for ve- 
locities and CaT EWs consists of 789 stars. We provide a 
detailed comparison between this data-set and the one from 
IWalker et ail l|2009l ) in the Appendix. 



2.1 Metallicity scale 

We use the CaT EWs to derive metallicities ([Fe/H]) for the 
target stars. 

The near-infrared CaT lines have been extensively used 
in the literature as empirical estimators of the [Fe/H] abun- 
dance of individual RGB stars observed at intermediate 
spectral resolution. 

This method has been empirically calibrated for 
stellar clusters and proven reliable over the range 

-2.1 < [Fe/H] < - .2 and 2.5 < a ge [Gyr] < 13 

(e.g. iRutledge et all Il997l ; ICole et al.l l2004h on the 
ICarretta fc Grattonl l| 19971 ) metallicity scale. Battaglia et alj 
1 2008?) tested the validity of the method for composite stel- 
lar populations, observing the same 129 individual stars 
in the Sculptor and Fornax dSphs both at intermediate 
and high resolution, and showed that the CaT EW- [Fe/H] 
relation can be applied to composite stellar populations 
over the explored rang e —2.5 < [Fe/H] < — 0.5. 

IStarkenburg et al, l|2010l ) carried out a synthetic spec- 
tral analysis of the CaT lines down to [Fe/H]= —4 in or- 
der to understand their theoretical behaviour as a function 
of physical parameters such as metallicity, gravity, effective 
temperature and possible biases in the low metallicity range. 
The authors find that the a simple linear relation no longer 
holds approximately below [Fe/H] < — 2.5. This behaviour 
is mostly a reflection of the change of the CaT lines from 
wing-dominated at high metallicity to core-dominated at low 
metallicity. 

The authors provide a revised relation, which holds over 
the much larger range — 4.0< [Fe/H] < — 0.5: 



[Fe/H] = -2.87 + 0.195 x (V - Vhb) + 0.458 x EW + 
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- 0.913 x Ely -1 ' 5 + 0.0155 x T.W x (V - Vhb) (1) 

In the above EW is the summed EW of the two strongest 
CaT lines, i.e. EW = EW2 + EW3, V is the apparent mag- 
nitude in V-band of the star and Vhb is the apparent mag- 
nitude in V-band of the horizontal branch of the galaxy. For 
Sextans we use Vhb = 20.35 from E rwin fc Hatzidimitrioul 
l|l995l ), which is also consistent with our photometry. 

Using our sample of 789 stars along the line-of-sight 
to Sextans , we find that t he [F e/H] values derived using 
Eq. (16) in Battagli a et al] (|2008l ) and Eq. fl} compare very 
well down to [Fe/H] = —2; at lower metallicities, Eq. (16) 
from iBattaglia et all (|2008h provides on average a larger 
[Fe/H] value than the revised calibration given by Eq. ([!}■ 
e.g. a [Fe/H]~ —2.4 with the old calibration corresponds 
to approximately [Fe/H]~ —2.6 with the revised one and 
[Fe/H] 2.6 to about [Fe/H] 3.0. 

In the following we will use Eq. fl]) for our CaT [Fe/H] 
determinations. In or der to derive the e r rors i n [Fe/H], we 
follow the approach o f lStarkenburg et al l (|2010^ l. i.e. we con- 
sider that the error due to the photometry is negligible with 
respect to the error in the determination of the CaT EW and 
we input in Eq. fl} the values of EW— Ioew and EW+Ioew 
and derive the corresponding la lower and upper limit in 
[Fe/H]. This reflects the fact that the relation is not linear 
and therefore symmetric errors in the EW determination can 
translate into asymmetric errors in [Fe/H]. 



3 MEMBERSHIP 

Given the low surface brightness of Sextans and its loca- 
tion on the sky, the RGB locus of Sextans is heavily con- 
taminated by Milky Way stars which are located along the 
line-of-sight to Sextans. In order to study the metallicity 
and kinematic properties of Sextans it is important that 
we first weed out these Galactic interlopers from our sam- 
ple. These interlopers will be mostly foreground MW stars. 
Given the Galactic coord inates of Sextans, the Besancon 
model ()Robin et al.ll200"3h predicts that the majority of in- 
terlopers along the line-of-sight to Sextans is made of dwarf 
stars, mostly from the thin and thick disk, and by smaller 
fraction by sub-giants and giant stars. 

Ideally one would like a method allowing a direct dis- 
crimination between members of the dSph (which have been 
selected to be RGB stars) from non-members (mostly dwarf 
stars) on the basis of, for example, spectral features indi- 
cating physical characteristics, such as in this case gravity. 
However, no such spectral features have been identified in 
previous studies using a similar spectral range and set-up. 
The common approach is therefore statistical. 

The simplest way of discerning between members and 
non-members is to consider the line-of-sight velocity distri- 
bution of the observed targets (see Figure[3]for the distribu- 
tion of heliocentric velocities for our VLT/FLAMES targets) 
and iteratively apply a k-cr cut around the systemic velocity 
of the galaxy, where a is the internal line-of-sight velocity 
dispersion of the galaxy. The k-cr cut in common use is the 
3a cut, justified by the close to Gaussian line-of-sight veloc- 
ity distribution observed for dSphs. 

We first follow this approach and then refine our selec- 
tion using further criteria. 



3.1 3-a clipping 

We first identify the peak of the velocity distribution and 
derive the weighted average velocity and dispersion for the 
stars with a velocity within 4a of the peak, using as a first 
guess a broad dispersion of a = 15 km s" 1 . We then re- 
peat the procedure restricting the range of selection to 3a. 
We then apply to the distribution of line-of-sight velocities 
for the stars withi n the 3a range the max imum likelihood 
method outlined in lHargreaves et~a l. (1994). This results in 
a systemic velocity w sys ,hoi = 226.3±0.6 km s _1 and a global 
dispersion Ohoi = 8.8±0.4 km s _1 for 182 probable members. 
The systemic velocity derived is in very good a greement with 
the va lue of ti sy s,hci = 225.8 ± 0.5 km s -1 from Walker et al] 
(2006); for the global dispersion they quote values from 3 
sample sizes: N= 276 <r hc i = 7.1±0.3 km s _1 ; N=294 <r hc i = 
8.9 ± 0.4 km s" 1 ; N=303 cr hc i = 10.3 ± 0.5 km s~ x . Our de- 
termination is in very good agreement with their determi- 
nation from the N=294 sample. We find a good agr e ement 
also with the values derived by lHargreaves et all (|l994l ) 
(«sys,hci = 224.4 ± 1.6 km s _1 , cr^el = 7.0t 1 ;f,km s -1 from 
21 members), and c onsistent values also to the study of 
ISuntzeff et all (|l993T l (v sys , hc i = 227.9 ± 1.8 km s _1 ; a heX = 

6.2 ± 0.9km s _1 based on 33 stars). 

Given the large extent on the sky of Sextans, the So- 
lar motion and Local Standard of Rest (LSR) will con- 
tribute a component along the line-of-sight to the object, 
which may introduce spurious velocity gradients and affect 
the derivation of quantities such as internal rotation for 
example. In order to remove this effect we transform the 
heliocentric line-of-sight velocities into line-of-sight veloci- 
ties in a frame at rest with respect to the Galactic centre 
(«gsr where GSR stands fo r Galactic standard of rest) . For 
this we use the formula in iBinnev fc Tremaind (|1987l ). us- 
ing a LSR velocity "Ulsr = 220 km s _1 at the Solar radius 
(R e = 8kpc) and a Solar motion of (U, V, W)=(10, 5.25, 
7.17) km s _1 (|Dehnen fc Binnevl[l99"sh , where U is radially 
inward, V positive in the direction of Galactic rotation and 
W towards the North Galactic Pole. 

We repeat the determination of the systemic velocity 
and global dispersion using the GSR velocities, and we ob- 
tain w S ys,GSR = 78.6 ± 0.7 km s _1 and a gsr. = 8.8 ± 0.5 
km s _1 for 182 probable members. Here after we will use the 
GSR velocities, except when comparing the individual line- 
of-sight velocities from this study those from other studies 
to ease the comparison. 

3.2 Use of the Mglline at 8806.8 A as a 
dwarf/giant discriminator 

Figure U shows the distribution of [Fe/H] measurements for 
targets along the line-of-sight to Sextans against their line- 
of-sight velocity in the GSR systerrQ. Around the systemic 
velocity of Sextans ± 3ctgsr one can distinguish the stars 
that are probable members of the galaxy, from those likely 

1 Note that the CaT relation relies on the assumption that all the 
stars can be considered at the same distance. While this is a rea- 
sonable assumption for the stars belonging to the Sextans dSph, 
this is most likely not the case for the foreground/background 
MW stars. Therefore, the [Fe/H] values for the stars which are 
probable non-members cannot be considered as meaningful. 
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to be foreground/background MW stars. From Fig. [4] it is 
clear that a simple kinematic selection like the 3-cr clipping 
is able to remove the majority of the MW contaminants but 
that a small fraction may still be present in the 3-cr region of 
membership selection. The fraction of contaminants appears 
to change with the [Fe/H] value derived from the CaT lines: 
at [Fe/H] < —2.0 no considerable contamination is expected, 
while at larger metallicities Sextans stars are more likely to 
be contaminated by Galactic stars whose velocity falls in the 
3cr membership selection region. 

Here we attempt a new empirical approach to refine our 
membership determination, using the Mglline at 8806. 8 A 
as a further possible discriminator between MW and Sextans 
members. 

We first select the stars that have a high probability of 
membership to Sextans by restricting the velocity range of 
selection to 2ctgsr just for the purpose of this analysis; the 
highly likely non-members are chosen to be the ones with 
GSR velocity more than 4<tgsr away from Sextans systemic 
velocity. Given the small EW of the Mg I line, in general less 
than 1 A , only spectra with S/N/A > 20 are used to explore 
possible trends without being affected considerably by the 
noise. In order to explore the behaviour of the Mglline with 
gravity in an empirical way, in Fig. [5] we show the behaviour 
of the EW of the Mg I line versus the CaT T.W in bins of V- 
magnitude, for both high likely members and non-members. 

As it is visible in Fig. [5] the loci of Mg I EW are quite 
distinct for members and non-members: members cluster 
at EW between ~0 and 0.2 A , while non-members occupy 
preferably the region of MglEW between 0.5 and 0.9 A. 
Some members have negative MglEW: this is just a re- 
flection of the lower S/N of those spectra; considering the 
error-bars these are consistent with MglEW close to zero, 
therefore with the line being hidden in the noise. 

It is clear that in general the non-members display a 
larger MglEW with respect to members. Preliminary use 
of synthetic stellar models suggest that this difference may 
not be driven exclusively by gravity but also of the different 
metallicity of the dSph star and MW stars ( with the former 
being below [Fe/H]< —1 and the latter peaking around —0.5 
to ~0). 

From the figure it appears that the trend of MglEW 
versus CaT TiW changes slightly among the bins in V- 
magnitude, which on a first approximation indicates a de- 
pendence on the gravity. However, a safe selection is to con- 
sider as members the stars with MglEW < 0.5 A . Hereafter 
we add this constrain to our selection criteria for probable 
Sextans stars. We have checked that this selection is valid 
also when exploring the behaviour of the MglEW versus 
CaT T,W in bins of V— I, as a proxy for effective tempera- 
ture. 

Applying this selection, 6 stars classified as members 
from the 3cr clipping appear to be contaminants (these are 
indicated as diamonds in Fig. 2} . 

We will improve on the simple empirical approach pre- 
sented in this section using synthetic spectral libraries in 
order to assess the applicability range of this method. 

3.3 Location on the CMD 

Our VLT/FLAMES targets have been chosen to have mag- 
nitudes and colours broadly consistent with that expected 



for RGB stars. We can now apply the membership criteria 
of the previous two sections to redefine the RGB locus oc- 
cupied by probable Sextans members. Figure [5] shows the 
location on the CMD of the probable members and prob- 
able non-members, selected using the criteria of Sects. 13.11 
and 13.21 All the probable members are closely located along 
the RGB ridge. The only evident outlier is found at VT~ 0.9 
and V~17.2 (see Fig. [6}. This star is located at a projected 
distance R ~2.2 deg. Given its metallicity of [Fe/H] = —1.9, 
its location on the CMD appears inconsistent with stellar 
evolution models; we therefore exclude this star from subse- 
quent analysis. 

3.4 Final sample of probable members 

By adopting all the membership criteria described above, 
i.e. 3cr velocity clipping, MglEW and location on the CMD, 
the final sample of probable members is 174 stars. We have 
also excluded one star having an error in [Fe/H] > 1 dex. 
Using this sample the systemic velocity and global dispersion 
in the GSR system are « sys ,GSR = 78.4 ± 0.6 km s and 
""GSR = 8.4 ± 0.4 km s" 1 , while in the heliocentric system 
they are « sys .hci = 226. 0± 0.6 km s _1 and a global dispersion 
o"hci = 8.4±0.4 km s _1 , all in good agreement with previous 
determinations. Hereafter we will use this sample for the 
analysis of the properties of Sextans. The location on the 
field-of-view of the cleaned up sample and the contaminants 
is shown in Fig. [7] From this plot the problems in collecting 
large samples of individual stars in an object like Sextans 
are clear: it is very extended, therefore many pointings of 
wide-field area spectrographs like FLAMES are needed, and 
because of its very low surface brightness, the outermost 3 
pointings did not yield any probable members. 



4 RESULTS 

4.1 Metallicity properties 

Metallicity measurements of large numbers of individual 
stars are important since they can shed light on the pro- 
cesses that drive galactic chemical enrichment and by mea- 
suring the metal abundance in stars over a wide age range 
and across the galaxy one can recover how the metal con- 
tent built up in time and how the enrichment proceeded 
throughout the galaxy. 

Figure [8] shows the metallicity distribution of proba- 
ble members of Sextans as derived from our VLT/FLAMES 
data (left: as a function of elliptical radius; right: overall dis- 
tribution). The stars in Sextans cover a wide range of [Fe/H] 
values, with the majority covering values from about —3.2 to 
— 1.4, but reaching down to —3.8 and up to —0.1; the average 
value of the distribution is [Fe/H] avg = —1.9, with a scatter 
of 0.6 dex, and the median value is [Fe/H] mc d = —2.3. 

The left panel of Fig. [8] clearly shows that the metallic- 
ity properties of this galaxy change with the projected ra- 
dius from the centre: the inner parts of the galaxy are more 
metal rich than the outer parts. This kind of behaviour is 
not unique to Sextans but has been detected also in other 
Milky Way satellites, on the basis of CaT data. However, 
here in Sextans the behaviour is more enhanced than in the 
other systems, with a step like behaviour of the metallicity 
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distribution going from the region at R<0.8 deg to the re- 
gion at R>0.8 deg: while the inner parts (R<0.8 deg) show 
a broad distribution, covering the whole range of [Fe/H] val- 
ues, the stars in the outer parts (R>0.8 deg) appear to 
have all [Fe/H] values < — 2.2 (except for 2 stars out of 
19). The current data appear to suggest a [Fe/H] thresh- 
old ([Fe/H]= —2.2) in the outer parts, but more data would 
be required to confirm that this result is not due to small 
number statistics. 

It is likely that the two "metal-rich" stars at R>0.8 
deg are contaminants: their [Fe/H] value is very discrepant 
with respect to the rest of the distribution, with metallicities 
of about 1 dex above the largest [Fe/H] of the rest of the 
stars. These stars may be RGB stars belonging to the stellar 
halo of our Galaxy. Similarly, at R~ 0.5 deg, one star is 
present with [Fe/H]~ —0.1, more than 1 dex larger than the 
maximum [Fe/H] displayed by all the other probable Sextans 
members anywhere in the galaxy. Even though it is possible 
that these 3 stars may belong to the stellar halo of the Milky 
Way, these objects full-fill all our membership criteria and 
we have no means to distinguish them in a direct way from 
genuine Sextans members, therefore we will keep them in 
the analysis. 

As shown in Fig.[9]the metallicity distribution of R>0.8 
deg is well approximated by a Gaussian with best-fitting 
peak position at [Fe/H]= — 2.73±0.06 and dispersion 0.24± 
0.05 dex (reduced x 2 = 0.6). The same figure shows that 
the region at R<0.8 deg is well represented by a sum of two 
Gaussians: one having the same shape parameters -modulo 
the amplitude - as in the outer parts, and the other Gaussian 
with best-fitting peak position at [Fe/H]= —2.04 ±0.03 and 
dispersion 0.25 ± 0.03 dex (reduced \ 2 = 0.7). At R<0.8 
deg these two Gaussians do overlap over a metallicity range 
and at [Fe/H]~ —2.3 the contribution of the stars belonging 
to the "metal-rich" component is that of the "metal-poor" 
starfl We will come back to this point in Sect. 14.331 

Figure [9] also shows the presence of 31 stars with CaT 
[Fe/H] < —2.8, whose presence was not revealed when us- 
ing the CaT calibration based on globular clusters; 11 of 
these stars have [Fe/H]< —3, i.e. they could be classi- 
fied as extremely metal poor stars. The presence of stars 
of these metallicities excludes the hypothesis that Sextans 
was formed from a pre-enriched medium of [Fe/H]> —3 as 
iHelmi et all l|2006l ) put forward. 

Seven Sextans stars with CaT [Fe/H] < —2.7 have been 
followed up at high spectral resolution bv lAoki et al.l (|2009h 
and Tafelmayer et al. (2010), showing that the CaT calibra- 
tion holds down to these low metallicities. 



4.2 Velocity gradients 

Early studies on dSphs, based on datasets of dozens of line- 
of-sight velocities in the central regions of these objects, had 
not yielded statistically significant detections of velocity gra- 
dients in these galaxies and therefore dSphs have commonly 
been regarded as exclusively pressure supported systems. 



2 Here the Gaussian fit to the metallicity distribution at different 
distances from the centre is only meant to highlight the differences 
in the peak metallicity value between the inner and outer parts. 



The gathering of much larger datasets, with hundreds line- 
of-sight velocities spread throughout the face of dSphs, led 
to the detection of velocity gradients in several of these ob- 
jects. 

Figure [10] shows a velocity field for Sextans obtained 
from the probable members of the galaxy. There are hints of 
the presence of a velocity gradient along the projected major 
axis of Sextans, but from this plot it is also evident that, 
given the small number of stars distributed over a large area, 
more details, such as velocity contours, cannot be reliably 
obtained. 

We analyse the velocity trends of the probable mem- 
bers along different directions across the face of the 
galaxy. Given the relatively sparse spatial sampling, we ex- 
plore only 4 directions: along the major and minor axes 
(P.A.= 56°, 146°, respectively) and along 2 intermediate 
axes (P.A.= 101°, 191°). We consider a slit of variable width 
along those axes and explore the trend of the line-of-sight 
velocity in the GSR frame as a function of the linear dis- 
tance from the centre of the galaxy for those stars within 
the "slit". We decided not to bin the data, as they are too 
sparse. 

Figure[TT]shows the results for a slit width of 0.3 deg. As 
fiducial functional form for the gradient we use a straight- 
line, with the intercept fixed to the systemic velocity, and 
we perform an unweighed fit to the data points referring to 
the probable members. In order to derive the error in the 
fitted slope, we create 10000 mock velocity datasets in which 
the stars are as many as the number of probable members 
and with the same positions, but the velocities are drawn 
from Gaussians centred on the measured velocities and with 
dispersion given by the velocity errors and the global disper- 
sion of Sextans added in quadrature. For each of these mock 
datasets we repeat the unweighed fit to the velocities along 
the various axes; the error in the slope is given by the la 
level in the distribution of fitted slopes. For this slit width 
the detected slopes are: 8.5 ± 3.0 km s" 1 deg -1 along the 
major axis, —0.8 ± 4.2 km s" 1 deg -1 along the minor axis, 
0.4 ± 5.3 km s" 1 deg" 1 along the axis at P.A.= 101° and 
7.5±3 o km s" 1 deg" 1 along the axis at P.A.= 191°. 

In order to assess the significance of such detections we 
repeat the above procedure reshuffling the measured veloc- 
ities 10000 times, where this should be equivalent to elimi- 
nating any coherent motion and therefore simulating a non- 
rotating dSph with the same overall dispersion as measured 
for Sextans. We then calculate how many times the mock 
datasets give amplitudes larger than the measured ones. For 
a slit width of 0.3 deg the percentages are 0.4%, 58.4%, 
46.9% and 0.9% along the major, minor, P.A.= 101° and 
P.A.= 191° axes, respectively! 3 ]. The only statistically sig- 
nificant detections appear to be the ones along the major 
axis and along the axis at P.A.= 191°. If these velocity gra- 
dients are due to rotation about the minor axis of Sextans, 
we would expect to detect similar velocity gradients along 
the P.A.= 101° and P.A.= 191° axes, given that they are 

3 We checked whether the gradient along the P.A.= 191° axis 
may be driven by the two points at distance larger than 0.9 deg, 
by repeating the determination removing these two points. The 
gradient is still present, 6.3 ± 3.7 km s -1 deg" 1 , although the 
significance slightly decreases (the percentage of mock dataset 
giving as large a gradient increases to 4.5%). 
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symmetric about the minor axis; however, while the spatial 
coverage along the P.A.= 191° is quite homogeneous and 
symmetric about the centre of the galaxy, this is not the 
case along the P.A.= 101° axis, along which we have data 
only in the N-W part of the galaxy; it is therefore unclear if 
the discrepancy in the detection is due to the difference in 
spatial coverage. 

We find that, when increasing the slit width, for exam- 
ple of 0.4 deg, the amplitude of the gradients slightly de- 
creases, which suggests that the kinematics along a certain 
axis starts to be contaminated by the inclusion of more stars 
moving differently because placed at different distances from 
the centre. In this case the percentage of the mock datasets 
giving amplitudes larger than the measured ones is larger 
than for a slit width of 0.3 deg, but the gradients along the 
projected major axis and the P.A.= 191° axis remain sta- 
tistically significant. 

Another concern in the derivation of the velocity gradi- 
ent is the possible inclusion in the sample of contaminants 
from the Milky Way, or vice versa the possible exclusion of 
Sextans stars based on the 3a cut-off. However, none of the 
two points appear to be an issue in the analysis here carried 
out as it is visible from Fig. [TT]where we also show the non- 
members falling in the considered slits: from these figures it 
is clear that the fraction of unidentified contaminants still 
present in the 3a velocity range is likely to be negligible, and 
also that we are unlikely to have missed probable members 
because of the 3a cut-off. 

As mentioned already, with the current spatial sampling 
it is difficult to extract a velocity field and to determine a 
rotation pattern. We adopt a simple empirical approach and 
assume that v TOt = fc*d m m, where k = 8.5 km s" 1 deg -1 and 
d m i n is the angular distance of each star from the minor axis 
along a direction parallel to the major axis (with d m i n posi- 
tive above the minor axis and negative below the minor axis, 
the formula gives receding and approaching velocities above 
and below the minor axis, respectively). To assess how good 
a representation of the data this is, we subtract the v Iot from 
the stars classified as probable members, and re-derive the 
velocity field and the trends along the various slits. There 
appears to be a residual pattern in the velocity field from 
the "rotation-subtracted" velocities; indeed, velocity gradi- 
ents are still present along the various axes, but they do not 
appear to be statistically significant. For a slit width of 0.3 
deg the detected slopes are: —0.1 ± 3.1 km s" 1 deg -1 along 
the major axis, —LSt^g km s" 1 deg -1 along the minor axis, 
7. it 5% km s" 1 deg" 1 along the axis at P.A.= 101° and 
1.6±3.1 km s" 1 deg" 1 along the axis at P.A.= 191° . The per- 
centages of mock datasets that give amplitudes larger than 
the measured ones are large, 48.4%, 66.4%, 9.0% and 30.6% 
along the major, minor, P.A.= 101° and P.A.= 191° axes, 
respectively, i.e. the gradients do not appear to be highly 
statistically significant, meaning that the rotation pattern 
we assume is a reasonable representation of the data. 

It appears that a weak velocity gradient is present in 
Sextans, but the direction and amplitude of this gradient 
is difficult to constrain. No te that our determin ation differs 
from the determination by I Walker et al.l (|2008T l. who find a 
gradient of -2.1 ± 0.8 km s" 1 deg" 1 along a PA - 120°. 
Given that the velocities between the two studies compare 
well (see Appendix), this difference is most likely due to 
the different coverage and size of the sample, as well as the 



methodology used. A larger statistical sample and coverage 
of the galaxy would help quantifying this better; however, 
given the low surface brightness and large extent on the sky 
of Sextans, this may be a challenging task. 

4.3 Observed velocity dispersion profile 

4-3.1 All members 

We derive the l.o.s. velocity dispersion profiles for Sextans 
by binning the velocities of stars at similar radii using the 
sample of probable members as determined in Sect. 13.41 and 
deriving the average velocity and dispersion per bin using 
the procedure described in Hargreaves et al. (1994). 

In deriving the l.o.s. velocity dispersion profile we 
checked for the effect of a number of issues. 

First, we check the effect that the number of stars per 
bin has on the derived l.o.s. velocity dispersion profile by 
fixing the number of stars per bin to a constant value of 
N=20,25,30 and 35 (except for the last bin, which is allowed 
to have a lower number of points). We find that, even though 
the shape and amplitude of the l.o.s. velocity dispersion pro- 
files are consistent within the error-bars among all the cases, 
the binning with N=35 stars is the most reliable one because 
it produces less oscillations in the trend of average veloci- 
ties per bin, and that these velocities are consistent with the 
systemic velocity. 

We then explore the effect of binning the data in el- 
liptical annuli (of constant ellipticity and position angle 
(e = 0.35 and P.A.= 56°) and circular annuli, using a con- 
stant number of N= 35 stars per bin. In fact, one uses stars 
with a non-spherical spatial distribution as kinematic trac- 
ers of an underlying distribution which may be spherical (as 
well as oblate, prolate or even triaxial). We also use a con- 
stant number of N= 35 stars per bin when exploring the 
effect of subtracting the detected velocity gradient onto the 
determination of the velocity dispersion profile (hereafter 
we refer to the profile derived from the individual velocities 
from which the velocity gradient has been subtracted as the 
"rotation subtracted" profile). These results are shown in 
Fig. (left): the l.o.s. velocity dispersion profiles derived 
from the elliptical and circular binning, when subtracting or 
not rotation, all agree with each others within the error-bars. 
In the following we adopt the rotation subtracted velocities 
and the elliptical binning. 

The effect of using bins of variable width which in- 
creases with projected radius, instead than using a constant 
number of stars per bin, is shown in the right-hand panel of 
Fig. 1121 We explored 3 different binnings, which are labelled 
as "Binningl" , "Binning2" and "Binning3" in the right-hand 
side panel of Fig. I12p . All these choices of binning are com- 
patible with the trend shown when keeping the number of 
probable members per bin constant. In the following, When 
deriving the mass profile of Sextans we will use the l.o.s. 
velocity dispersion profile from "Binning2" ; however, the ef- 
fects of a different binning on the mass determination will 
also be considered, by checking the results obtained when 
using the profile from "Binning3". We will not use "Bin- 
ningl" for the analysis because it represents the case where 
the average velocity per bin differs the most from the sys- 
temic. 

The l.o.s. velocity dispersion profile of Sextans appears 
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to slightly increase with radius although, given the large size 
of the error-bars in the last bins, it is also compatible with 
remaining constant. 

To c ompare the l.o.s. ve locity dispersion profile pub- 
lished bv lWalker etafl (|2007f ) and ours we should consider 
the case of the circular binning, with no velocity gradient 
subtracted: in this case our profile remains approximately 
consta nt around 8 km s -1 , consistent with the determina- 
tion bv lWalker etafl (|2007f ). 

4.3.2 "Metal-rich" and "Metal-poor" stars 

In previous work, we have combined the information on the 
spatial location, metallicity and velocity of individual stars 
in dSphs. This allowed us to explore links between the metal- 
licity and the kinematics of stars in the Sculptor and Fornax 
dSphs (T04 and B06, respectively), and shown that in these 
systems metal-rich stars - centrally concentrated - have a 
smaller velocity dispersion than metal-poor stars - which in 
turn have a more extended spatial distribution. The metal- 
licity chosen to illustrate this different kinematic behaviour 
between the colder and hotter stellar components is some- 
what subjective, but it is chosen on the basis of consider- 
ations on the metallicity distribution in order to minimize 
contamination between the two metallicity components. 

As seen in Sect. 14. fl at R<0.8 deg, the region where the 
whole range of metallicities is present, the metallicity distri- 
bution of Sextans is well-approximated with the sum of two 
Gaussians, a "metal-rich" and a "metal-poor" one, which 
cross at [Fe/H]~ —2.3. Therefore to minimize contamina- 
tion between these two components, we exclude the range 
—2.4 < [Fe/H] < —2.2 and we derive the trend of mean 
velocity and dispersion for the stars more metal-rich than 
[Fc/H] = —2.2 and more metal-poor than [Fe/H] = —2.4 
over different spatial bins. 

Figure [T3] shows that also in Sextans "metal-rich" and 
"metal-poor" stars display a different behaviour (akin to our 
findings in the Fornax and Sculptor dSphs), with the former 
population in general kinematically colder than the latter. 
The profile of the metal-rich stars appears to be slightly de- 
clining, while for the metal-poor stars it is everywhere con- 
stant (except in the first bin, see next section) and extends 
further than the profile for the MR stars, as a consequence of 
the larger spatial extent of the metal-poor stars. It remains 
to be proven if the different spatial extent of the stellar pop- 
ulation analysed is the only factor that causes the different 
kinematics, or if one needs to take into account also differ- 
ent orbital properties for the different stellar populations. 
Also, the outer bin for the "metal-rich" stars suffers from 
low number statistics and would certainly benefit from a 
larger sample of stars with [Fe/H]> —2.2. 

4.4 Substructures 

iKlevna et afl (|2004r ) detected a kinematically cold substruc- 
ture, with dispersion close to zero, from 7 stars with radial 
velocity measurements in the inner 5 arcmin of Sextans. One 
of the hypothesis that they put forward is that the stars 
giving raise to this kinematic feature may belong to the 
remnant of a disrupted stellar cluster which spiraled into 
Sextans centre. The authors note that this hypothesis can 



account for the sharp central rise in the light distribution of 
Sextans (e.g. Irwin & Hatzidimitriou f995). Also the distri- 
bution of blue stragglers in Sextans, with the brighter (i.e. 
more massive) blue straggler s more ce ntrally concentrated 
than the fainter ones (|Lee et alj [2003D could be explained 
by this mechanism. The authors estimate a luminosity of 
1.3xl0 5 L Q for the disrupted cluster. The existence of this 
central cold substruc ture has not been co nfirmed by subse- 
quent studies such as I Walker et al.l (|2006t ). who instead find 
another cold substructure of estimated luminosity ~ I0 4 
Lq at a different location, i.e. around the core radius of the 
galaxy. 

The inner point in Fig. Q3] shows peculiar characteris- 
tics for the metal-poor stars: the velocity dispersion is very 
cold, f .4±f.2 km s _1 , and the average velocity is 72.5±1.3 
km s _1 , about 4a away from the systemicQ. In the binning 
used these values are derived from the 6 innermost metal- 
poor stars, however from panel (a) of Fig. [13] it can be seen 
that other 3 MP stars at R< 0.22 deg share very similar 
velocities to this group of stars; only I MP star at R< 0.22 
deg is found at very different velocities and can therefore 
be thought of as belonging to the main Sextans population. 
This group of 9 stars do not only share very similar distances 
and kinematics but also metallicities, with an average metal- 
licity of [Fe/H]= —2.6 and a scatter of 0.15 dex. Note that 
the average metallicity error for the stars in the substruc- 
ture is 0.29 dex, much larger than the measured scatter in 
the metallicity, which means that the intrinsic scatter is go- 
ing to be much lower than 0.15 dex. The similar metallicity 
of the stars found in the cold kinematic substructure would 
point to these stars previously belonging to a single stellar 
population, i.e. a stellar cluster. 

Note that the information on metallicity is useful for 
such kind of detections: when considering all the stars to- 
gether in the inner bin, with no distinction made on metal- 
licity, the velocity dispersion is much larger, consistent with 
the rest of the galaxy. 

Assuming that the ratio of stars in the substructure (9) 
with respect to the total number of probable Sextans mem- 
bers (174) is representative, then the substructure would 
account for 5% of the overall Sextans population. Assum- 
ing a luminosity for Sextans from Table [T] and a stellar 
mass -to-light ratio (M/L)i um ,y = 2 typica l of globular clus- 
ters (|lllingworthl 11979 ; iPrvor et al. I ll988T ) as plausible for 
an old stellar population, then crude estimates of the to- 
tal luminosity and mass of this structure are 2.2 x 10 4 
L©and 4. 4 x 10 4 Mp. I f this substructure is the same one 
found in IKlevna et~afl (|2004T ). our estimated luminosity is 
approximatel y one order of magnitude smaller than what 
IKlevna et afl i|2004h estimated, that is 1.3 x 10 5 L . How- 
ever the estimate w as carried out in a different way, where 
IKlevna et alj (|2004l ) assumed that all the light in the central 
10 arcmin is due to cluster debris. 

We can instead assume that the estimated light in the 
central 10 arcmin, 1.3 x 10 5 Lq, is actually given by the 
sum of the light from the substructure plus the light from 
the overall Sextans population. W e then use the value of 
the central surface brightness from llrwin fc Hatzidimitrioul 



Here with a we denote the error in the determined average 
velocity for the substructure. 
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1 19951 ) - derived excluding the central point giving raise to 
the sharp central rise in the light distribution of Sextans - to 
estimate the contribution of the overall Sextans population 
to the light in the central 10 arcmin: this is about 7.8 x 10 4 
Lq. Therefore the contribution of the substructure would 
then amount to 5.2 x 10 4 Lq, much closer to our estimate. 
It is therefore possible that the s ubstructure here de tected 
is the same one as in the work of IKlevna et al l \2004 ). 

The properties we find are compatible with other cold 
substru ctures found i n oth er dSphs, such as in Ursa Mi - 
nor by IKlevna et all (|2003h . in Sculptor (|Battaglial I2007T ) 
and at projected distanc es of about 1 core radius in Sex- 
tans ll Walker et al.l [20061. This luminosity corresponds to 
an absolute My = —6.2, which is similar to globular clus- 
ters found in the Fnx, Sagittarius dSp hs and in other dwarf 
galaxies (e.g., see Ivan den B orah 20o|). 

To date, the most metal-poor known globular cluster 
resides in the Fnx dSph, w ith high resolution m easurements 
yielding a [Fe/H] = -2.5 l|Letarte et all 120061 ). If the stars 
in the cold substructure here detected do indeed belong to 
the remnant of a disrupted globular cluster, the [Fe/H] val- 
ues we measure would place it among the most metal poor 
globular clusters known. Given the large error-bars of the in- 
dividual CaT measurements for these stars, determinations 
from high resolution spectroscopic measurements would be 
needed to place this on a secure foot. Detection of the O-Na 
anti-correlation would confirm the hypothesis. 

Note that the CaT [Fe/H] values are calculated assum- 
ing that the stars in the substructure are at the same dis- 
tance of Sextans. If those stars are actually an external fea- 
ture, then the [Fe/H] values would not be correct. However, 
given the average velocity of the substructure well within 
1o"gsr from the systemic velocity of Sextans, it is likely that 
these stars do indeed belong to the dSph. 



5 MASS DETERMINATION 

We determine the mass of Sextans by means of Jeans mod- 
eling of a spherical and stationary system assuming different 
dark matter density profiles. 

The methodology we use consists of comparing the ob- 
served l.o.s. velocity dispersion <7i os for each distance bin 
with that predicted from the Jeans modeling for different 
models of the dark matter density profile and hypothesis on 
the velocity anisotropy of the traceiQ (see Sect. 15. 1]) . We ex- 
plore the space of parameters which define each model and 
determine the \ 2 as: 



2 

X = 



E 



Closi — CTlos(Ri;Pl3, 



(2) 



The variable p denotes a characteristic mass parameter of 
each dark matter model, pp denotes a parameter describ- 
ing the behaviour of the velocity anisotropy (see Sect. 15.2]) . 
Finally, 6j is the error in the observed l.o.s. velocity disper- 
sion. The best-fitting parameters are defined as those for 

5 By "tracer population" we mean those objects whose kinemat- 
ics can be used to recover properties of the total potential. In 
this case our tracers are the spectroscopically observed red giant 
branch stars. 



which x is minimized. We quote as errors in the individual 
parameters the projections of the A\ 2 = 2.3 region (corre- 
sponding to the region of 68.3% joint probability for a two 
free parameters \ 2 distribution). 

5.1 Predicted velocity dispersion profile 

The Jeans equation for a stationary spherical system in ab- 
sence of net streaming mo tions in any of the directions is 
l|Binnev fc Tremaine|[l987l ): 



1 d(p*ar,*) + 2/3cr; 



dr 



d4 
dr 



v 2 



(3) 



where p* is the density of the tracer; f3 is the velocity 
anisotropy parameter, defined as /3 — 1 — cr| „/crj? », assum- 
ing erg „ = °V,*> an d cr</>,» are the velocity disper- 
sions in the (r, 9, <p) direction respectively; (j> an d V c are the 
potential and the circular velocity of the total mass distri- 
bution. Note that /3 = if the velocity ellipsoid is isotropic, 
P = 1 if the ellipsoid is completely aligned with the radial 
direction, and /? < for tangentially anisotropic ellipsoids. 

The quantity to compare to the observations is 
the l.o.s. velocity disp ersion of the tracer population 
l|Binnev fc Mamonlfl982T ): 



Clos(-R) 



P*( r )< j2 ,* r 



R 2 



(l-p^)dr 



(4) 



where R is the projected radius (on the sky) and E*(ii) is 
the mass surfac e density of the t r acer. 

We refer to lMamon fc Lok as (2005) for the derivation of 
the l.o.s. velocity dispersion profile from the Jeans equation 
using different hypotheses on j3. 

As the discovery of multiple stellar populations in dSphs 
is a recent one, traditionally dSph galaxies have been tre ated 
as single component systems. In iBattaglia et alj (|2008T ) we 
modelled each stellar component separately in the Sculptor 
dSph, showing that this allows to relieve some of the degen- 
eracies. This kind of analysis needs imaging such that the 
surface number density profiles of the different components 
can be accurately determined. This is not the case to date for 
Sextans, therefore until better photometric data will come 
we will treat all the stars in Sextans as in one component. 

Below we discuss in detail the ingredients to solve the 
Jeans equation. 

• The spatial distr ibution of the tracer: 
llrwin fc Hatzidimitrioul l|l995h found that the surface 
brightness profile of Sextans is best fit by an exponential 
profile with scale radius R e — 15.5'. This is equivalent to 
assuming a Sersic profile with Sers ic radius Rs = 15.5' and 
shape param eters m = 1 as done bv lWalker et al.1 (|2007f ) and 
Lokas (2009). For the mass-to-light ratio we use (M/L)y = 1 
(M/L)v;0 and luminosity and distance from Table [1] 

The 3D density profile is derived from the surface 
brightness profile through inversion of Abel integrals, as- 
su ming that the star s are spherically distributed. We refer 
to lLokas et all (|2005h for all the relevant formulas. 

• The kinematics of the tracer - Velocity gra- 
dient: The expression for the projected velocity dispersion 
(Eq. [4| , which measures the projected random motion in a 
galaxy, has been derived for the hypothesis that the system 
is not rotating. In systems where the rotation law can be ac- 
curately derived and there is enough statistics to derive the 
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dispersion profile along both the major and minor axes, one 
can use the Jeans equations in cylindrical coordinates. The 
present coverage of this dataset does not allow this alterna- 
tive. Therefore we use the observed l.o.s. velocity dispersion 
profile derived from the rotation subtracted individual ve- 
locities. 

• The kinematics of the tracer - Velocity 
anisotropy: Since the variation of the velocity anisotropy 
with radius is not known (as this requires proper motions), 
we consider two hypotheses: /3 constant with ra dius; and 
using the Osipkov- Merritt parametrization for ft (|Osipkovl 
ll979l : lMerrittlll985h 



/3om = r /(r +r a ) 



(•») 



where r a is the anisotropy radius. In the Osipkov-Merritt 
parametrization the anisotropy is always > 0, i.e. it is never 
tangential. The central regions are isotropic, and for r a — > oo 
the anisotropy becomes purely radial. At r — r a , P — 0.5. 
The smaller r a , the faster the anisotropy becomes very ra- 
dial. Models with large r a correspond to models with almost 
isotropic behaviour. 

• Total mass distribution: Since from previous works 
the contribution of the stars to the overall potential has been 
shown to be negligible in Sextans and in dSphs in general, 
we consider the dark matter halo as the only contributor to 
the kinematics of the tracer population. 

We consider two different models for a spherically sym- 
metric dark-matter halo potential: 

- Pseudo-Isothermal sphere. This model has been ex- 
tensively used in the context o f extragalactic rotation curve 
work (see S wa ters et al J 12000. and references therein). The 
density profile associated to this model is: 



P(r) = Po 



where r c is the core radius, po — 



K 2 (oo) 



(6) 



is the central den- 



AnGrl 

sity and V c (oo) is the asymptotic circular velocity. At large 
radii the density behaves as p oc r~ 2 . 
The resulting mass profile is: 



M(< r) = 47rp r 2 (r 



■ arctan — 
rJ 



(T) 



The profile is completely defined by po and r c , or any couple 
of non-degenerate parameters. 

- NFW model. This profile is motivated b y cosmological 
N-body sim ulations in a CDM framework ( Navarro et al.l 
1 19961 . 1 19971 ). In this case the DM density profile is given by 



p(r) = 



5 c p° 



(r/r„)(l + r/r s 



(8) 



where r B is a scale radius, p° the present critical den- 
sity and S c a characteristic over-density. The latter is de 
100 c 3 g{c) 



fined by 8 C — 



where 



= r v /r s is the con- 
its virial radius, and 



centration parameter of the halo, r- 

g(c) = (ln(l + c) — c/(l + c))~ 1 . The concentration c has been 
found to correlate with the halo virial mass M v in the range 
10 11 -10 14 ft~ 1 MmllNavarro et alj Il997l ; iBullock et all l200ll : 
IWechsler et alj|2002h . so that at a given redshift more mas- 
sive haloes have lower concentrations. In principle, the rela- 
tion between c and M v makes the NFW model completely 
defined by one parameter (e.g., M v ). However at a fixed 



mass the scatter in the predic ted concentration is large, of 
the order of A(log 10 c) = 0.18 (|Bullock et alj|200ll ). thus we 
do not consider the NFW density profile as a one-parameter 
family, but we describe it both by the concentration c and 
by the virial mass or the circular velocity at the virial radius. 

At large radii (for r 3> r s ), the density behaves as p oc 
r -3 , and therefore, the total mass diverges logarithmically. 
The resulting mass profile is 



M{< r) = M v 



1 + x 



(9) 

l|Klypin et alj 



where x — r/r s and f(x) = ln(l + x) 
12002ft , 

When integrating Eq. [4] we set the upper integration 
limit to r v where we use r 2/9 p*a"^*| rv = (we are essentially 
assuming that the particles are bound out to the virial ra- 
dius). As we will see in the next Section and as it is found in 
previous works on dSphs, the extent of the luminous matter 
is typically one order of magnitude smaller than the extent of 
the best-fitting DM halo, thus setting the upper integration 
limit to r v instead of to infinity will not affect our results, 
only decrease the computation time. 



5.2 Results 

We perform the fitting procedure for the observed l.o.s. pro- 
file derived in Sect. l4~3TTl i.e. "Binning2" 0. The parameter 
p in Eq. denotes the mass enclosed within the last mea- 
sured point for the isothermal sphere and the virial mass for 
the NFW model; the parameter pp is pg = /5 for the models 
using /3 constant with radius and pp = r a for the models 
using Pom- 

The results are summarized in Tab. [3] 

• Isothermal sphere: In this model the total mass of 
the DM halo is not finite, therefore we choose to fit the value 
of the mass within the last measured point, M(< i?i as t), 
which for "Binning2" is i?i aa t =2.36 kpc. The other free pa- 
rameter to take into account is the velocity anisotropy p. 

We explore the performance of models with core radii r c — 
0.001, 0.05, 0.1, 0.5, 1, 1.5, 2, 2.5 and 3 kpc. In practise, 
we fix the core radius to each of these values and mini- 
mize the x 2 to obtain the best-fitting M(< i?i as t) an d /3. 
The results of the fits are summarized in the plots shown 
in Fig. [14] (left). At a given core radius, the best fit mod- 
els with constant /3 give in general smaller \ 2 values than 
when P = Pom- The models for P — Pom and r c < 0.5 
kpc have a reduced \ 2 > 3, which for 4 degrees of freedom 
corresponds to a significance level of less than 2%; they can 
therefore be considered as statistically unacceptable. All the 
other best fitting models for the explored core radii have re- 
duced x 2 ^ 1) therefore they are all statistically acceptable. 
For these models the best- fitting M(< i?i as t) ranges from ~1 
- 4 xlO 8 Mq (smaller core radii, smaller masses); the best- 
fitting anisotropics ranges from —6 to 0.5 for P = const 



6 We checked that using the "Binning3" velocity dispersion pro- 
files produces best-fitting parameters and models that are very 
similar to those obtained when adopting the "Binning2" profile, 
although in this case the minimum x 2 values for the best-fitting 
models are larger than the ones from "Binning2". 
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(smaller core radii more tangential anisotropies) , while the 
anisotropy radius goes from 2 to 15 kpc (the upper limit of 
the explored r a ), forcing the anisotropy to be close to zero. 
The models with the smallest reduced x 2 (see Fig |15[) for 
P = const have r c = 3 kpc (M(< R laa t) = 4 ± 0.8 x 10 s 
M and 0.06 < P < 0.6) and for P = Pom r c = 1.5 kpc 
(M(< ftast) = 3.2 ± 0.7 x 10 8 M Q and r a >1.5 kpc). 

This tendency for increasing tangential anisotropy with 
decreasing core radius comes from the fact that smaller core 
radii models have higher concentrations of mass in the cen- 
tral regions and this tends to increase the central value of 
the velocity dispersion. To match a nearly constant velocity 
dispersion profile the anisotropy then needs to be tangential. 

• NFW model: For this model we let the concentration 
vary from c =10, 15, 20, 25, 30, 35, and for each of these 
values we derive the best-fitting virial mass and velocity 
anisotropy. Using an e xtrapolation of t he formulae in th e 
N-body simulations of Ijind (|2000l ) (see iKoch etafl l2007h . 
the concentrations adopted here would correspond to the 
concentrations expected for haloes with virial masses in the 
range 10 6 -10 12 M Q (smaller mass larger c) well containing 
the range of virial masses expected for dSphs. 

At fixed concentration, all the models with /9 = Pom have 
reduced x 2 larger than then models with P = const. The 
models with c > 20 have a reduced \ 2 > 2.5, correspond- 
ing to a significance level of less than 5% for 4 degrees of 
freedom and have therefore low statistical significance. The 
best result for the P = const case is also given by the model 
with c = 10, with a virial mass M v = 2.6±0.8x 10 9 M and 
P = -0-GtiA ( see FigHS). All the models with P = const 
have reduced \ 2 between 1 and 1.5, where the virial masses 
range from approximately 3x10 s to 2.6xl0 9 Mq (smaller c 
larger virial mass) , which gives a mass within the last mea- 
sured point of l-2xl0 8 Mq, and P goes from about —4.5 
to —0.6 (smaller c less tangential p). This also explains why 
most of the models with P — Pom are not favoured for NFW 
profiles: in the hypothesis of P — const, which yields a good 
fit to the data, the best-fitting anisotropies are mildly tan- 
gential; since in the P = Pom the anisotropy is forced to 
be positive, than the resulting model cannot reproduce the 
data as well. 

We conclude that because of the well-known mass- 
anisotropy degeneracy, it is not possible to firmly distinguish 
between a cored and a cuspy dark matter distribution, nor 
about different hypothesis of anisotropy, solely on the basis 
of the line-of-sight velocity dispersion profile. The only mod- 
els that can be excluded among the explored ones because 
of their large \ 2 values are the NFW models with p = Pom 
and c > 20. However, the general trend is for the best-fitting 
models to have P close to isotropic, be it slightly tangential 
as for the NFW case or slightly radial as for the cored profile. 
Also, the tendency for cuspy profiles of low concentrations to 
produce better fits than for larger concentrations seems to 
suggest that shallower profiles are preferred. Independently 
on the model, we find that the mass within the last mea- 
sured point of Sextans, at 2.3 kpc, should be in the range 
2-4 x 10 s Mq. This would give a dynamical mass-to-light ra- 
tio between 460 and 920 (M/L)y )0 . 

The presence of the central substructure does not affect 
the mass estimate derived. Although the velocity dispersion 
for stars at the distance covered by the substructure, 7.6±1.1 



km s , is lower than if these stars had not been included, 
i.e. 8.6±1.4 km s _1 , these two values are consistent within 
the errors. 



6 DISCUSSION 

• Metallicity properties: The stars in Sextans show a 
wide [Fe/H] range, with the majority having [Fe/H] values 
between —3.2 and —1.4, with an average [Fe/H] = —1.9. 
Among the sample of probable members, 11 extremely metal 
poor stars are found, reaching down to [Fe/H] = —3.8. If the 
stars in Sextans were formed from a pre-enriched medium, 
then the presence of stars with [Fe/H] as low as —3.8 sug- 
gests that this interstellar medium most likely had not been 
pre-enriched to a larger value. 

Sextans exhibits a clear spatial variation of its metal- 
licity properties: the region at projected radius R < 0.8 
deg displays the whole range of metallicities, while at larger 
projected radii practically only stars more metal poor than 
[Fe/H]~ —2.2 are present. Photometric studies of its stellar 
population show an age difference of at least 3 Gyr between 
the central and outer parts of Sextans, with the central parts 
containing stars between 10 and 14 Gyr old and the outer 
parts mostly 14 Gyr old stars. It appears then that the outer 
parts were formed very quickly, reaching a metallicity up to 
[Fe/H]— 2.2, while the inner parts could continue their evo- 
lution for at least 3 more Gyr and enrich to [Fe/H]~ —1.4. 
One can speculate that the less prolongued star formation of 
the outer parts may have been caused either from external 
processes, which removed the gas from the external region, 
or to internal effects, which depleted the outer parts of gas. 
It is although possible that gas was still present in the outer 
regions for some time but simply that conditions for star 
formation were not met. 

Similar, although less pronounced, spatial variations of 
metallicity properties have been detected in other d Sphs, 
such as for example in Sc ulptor |Tolstov et al.| [2004). For- 
nax dBattaglia et al.ll2006t ) and to a minor extent in Carina 
(|Koch et al.ll2006|^ There are hints that Leo I may exhibit 
such variations (jciullieuszik et al.l I2001S1 ) although the cur- 
rent data are not sufficient to put thi s result on a firm 
ground. Using Stromgren photometry iFaria et al.l (120071 ) 
show the presence of a possible spatial variation in the metal- 
licity properties of the Draco dSph. It is still unclear what is 
the driving process behind these metallicity variations. To 
disentangle environmental effects from internal evolution it 
is important to obtain observations of the stellar popula- 
tions, kinematics and metallicity variations of isolated dwarf 
galaxies, outside the halo of large spirals. 

• Velocity gradient: A statistically significant veloc- 
ity gradient is found along the projected major axis and an 
axis at P.A .= 191° in our samp le of 174 probable members 
of Sextans. Walker et al.l l|2008h also detect a velocity gradi- 
ent in Sextans, although with lower amplitude and approx- 
imately along the minor axis of the galaxy. Pinning down 
the amplitude and direction of such gradients, expected to 
be small for dSphs, would require large sample of stars not 
only in the central regions of the object but also in the outer 
parts, so that the trends along various directions may be ac- 
curately determined. Unfortunately this is very challenging 
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for a system like Sextans, because of its large extent on the 
sky and the outer parts being sparsely populated. 

Also a univocal interpretation of this feature does not 
appear possible at this stage. One possibility is that the gra- 
dient may have originated as the result of tidal disruption of 
Sextans from the Milky W ay, as predicted in N-body models 
of disrupted satellites fe.g. lOh et al.|[l995h . However no clear 
signs of tidal disruption such as tidal tails or S-shaped con- 
tours are found in Sextans that could confirm this hypothe- 
sis, although this may be just due to t he observational di f- 
ficulties in detecting such features (e.g. iMunoz et al] [2008). 
especially in a diffuse and heavily contaminated object like 
Sextans. Knowledge of the accurate proper motion of the 
galaxy would allow to reconstruct its orbital history and de- 
termine if strong interactions with the Milky Way may have 
occurred. However no direct determinations of the proper 
motion of Sextans are available yet in the literature. Such 
determinations would also allow to understand whether the 
observed gradient may merely be a geometrical effect due 
to the projection of the object transverse motion along the 
line-of-sight because of its large extent on the sky, since in 
this case the direction of the maximum gradient is expected 
to be aligned with the proper motion direction. 

Another possible option is that the velocity gradient 
is due to intrinsic rotation of the object. In this case, if 
the flattening of Sextans is due to intrinsic rotation, then it 
would be reasonable to expect the maximum gradient along 
the projected major axis of the object, as detected in our 
data. When rotation is responsible for the flattening of an 
axisymmetric galaxy, a simple relation holds between the 
true ellipticity of the galaxy and its mass weighted rotational 
velocity over the mass we ighted velocity dispersion, vo/ao 
jBinnev fc Tremamdll987h . As for estimate of the rotational 
velocity we use the fitted value of the gradient at R ~ 0.7 
deg, the distance of the outermost stars along the major 
axis which cause the rotation signal in our data-set: this 
value is w ro t,ios = 5.9 km s _1 and it is likely to represent a 
lower limit to the maximum rotational velocity. Replacing <7o 
by the velocity dispersion measured along the line-of-sight, 
and n by Wrot.los, the ratio « ro t,ios/o~GSR = 0.7 is very close 
to the value expec ted for an "isotropic rotat or" (0.727) of 
ellipticity e =0.35 l|Binnev fc Merrifieldl[l998l . as tabulated 
in the). This means that the shape of Sextans is consistent 
with being flattened by rotation, although due to the fact 
that Sextans velocity ellipsoid and inclination are unknown, 
it is possible that velocity anisotropy is also contributing to 
the flattening of the galaxy. 

At the moment it is still difficult to build a coherent 
picture not only on the cause of velocity gradients in Milky 
Way dSphs, but of eve n their presence. No significant gradi- 
ent is found in Leo II (iKoch et al.ll2007ft where the coverage 
is very homogeneous. [Battaglia et" al. (2008) found a signifi- 
cant velocity gradient along the projected major axis of the 
Sculptor dSph, confirmed by [Walker et al.l (|2009T l . Gradients 
of smaller amplitudes have als o been detected in Carina and 
in Fornax (Wa lker et al-l feoOO). in this case along directions 
close to the projected minor axis of the object. For those two 
objects the amplitude and direction of the gradient are con- 
sistent with being due to the transverse motion of the object 
along the line-of-sight, given their measured proper motion. 
However, while Carina has a very homogeneous coverage, 
this is not the case for Fornax, where the outer parts still 



await careful scrutiny. The outer parts of these objects ap- 
pear to be particularly important in the measu rement of ve- 
locity gradients, as shown by the case of Leo L lMateo et al.l 
(2008) found a highly statistically significant gradient when 
considering only the stars located at projected distances 
larger than 400", while they found no significant gradient 
neither when using the sample of stars at smaller projected 
distances, nor when considering the whole sample (which is 
dominated in number by the stars at R < 400"). Also com- 
bining data-sets together may give different resul ts: when 
using only their CMOS dataset lKoch et all (|2007l l find no 
statistically significant gradient in Leo I, while when combin- 
ing their CMOS dataset to their DEIMOS one the statistical 
significance of the detected gradient increases. 

This wealth of results on the different Milky Way dSphs 
shows that a coherent picture of velocity gradients as in- 
trinsic feature in dSphs is still lacking, and so probably are 
datasets with both large enough statistics and homogeneous 
coverage to pin down values and directions of these gradi- 
ents. 

Velocity gradients have been detected along the pro- 
jected major axis o f the 2 isolated dSphs of the Local Group : 
Cetus and Tucana dLewis et al.1120071 : iFraternali et alj|2009 . 
respectively j3- At the distance of Cetus and Tucana, un- 
less o f invoking three-body interactions or such (ISales et al.l 
|2007|) . tidal stripping from either the Milky Way or M31 
cannot have been important, therefore it is very unlikely 
that the detected gradients may be attributed to tides. Also, 
gradients due to the transverse motion of the system are ex- 
pected to be negligible due to the small extent on the sky 
of these objects. This would point to velocity gradients as 
being an intrinsic feature of these systems. 

Also other spheroidal systems in the Local Group do 
show velocity gra dients along their projected major axis. 
iGeha et all (|201Cj ) found significant rotation in the NGC 147 
and NGC 185 dEs. The authors calculate that the flattening 
of these galaxies is consistent with the amount of rotation 
found, although part of it may be due to anisotropy. These 
galaxies are dEs, not dSphs, with luminosities of the order 
of 10 s Lq, about 2 orders of magnitude larger than the 
common dSph, and 1 order of magnitude large r than For- 
nax. However , as d iscussed in iKormendv et al l (2009) and 
iTolstov et al] (|2009l ), the dEs in the Local Group (NGC 147, 
NGC 185 and NGC 205) share similarities with the dSphs, 
suggesting that these may be the high luminosity end of the 
sequence. 

• Mass determination: We have performed a Jeans 
analysis of the observed line-of-sight velocity dispersion pro- 
file of Sextans and find that for the best fitting cored profiles 
the mass within the last measured point ranges from 3.2- 
4xl0 8 Mq, while the mass within the last measured point 
for the best fitting NFW profile is smaller, 1.9 ± 0.6 x 10 s 
M©but consistent within la with the determination from 
the cored profiles. This mass range implies that Sextans con- 
tains between about 450-900 times more dark than luminous 
matter, in the typical range of Milky Way dSphs. 

Other authors have estimated the mass of Sextans 
at smaller distances from the centre. The mass that 



7 Note though that in these two cases directions other than the 
major axis have not been well explored. 
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IWalker et all ([20071 ) find within their last measured point 
at 1.1 kpc is 5.4xl0 7 M for a constant anisotropy and an 
NFW halo (with best-fitting M v = 3 x 10 s M Q and c = 20 
deriv ed from the concentration-virial mass relation in Ijind 
|2000| and best-fitting /3 = —2). This is consistent with our 
model with c= 20, which yields M v = 6±?"g x 10 8 M Q and 
P = -1.7 and M(<1.1 kpc)= 5 5±? j x 10 7 M Q . Note that 
in their analysis IWalker et al.l l|2007h choose to fix the con- 
centration to the virial mass according to the formula of 
IJind |2000l ) instead of exploring various concentrations for 
a given virial mass like in this work. It is possible that also 
IWalker et al.l |2007l ) data would have yielded a lower best- 
fitting concentration, if the correlation M v -c was not as- 
sumed. 

It has been recently shown that the determination of 
the dark matter mass wi thin a certain distan ce is robust to 
model assumptions (e.g. IStrigari et alj [20071 '). An interest- 
ing result has been that the mass enc losed within 0.3 (0.6 
kpc) spans a narrow ran ge of values l|Gilmore et al.l [20071 ; 
IStrigari et al.l [20071 . I2008T ) and this could imply a minimum 
mass scale fo r the existence of luminous satellites (see also 
lMatedll998l. for a sim ilar suggestion). In the determination 
of lStrigari et all |2007l ) the mass enclosed within 0.6 kpc for 
Sextans is M(<0.6 kpc)= 0.9±o:3 x 10 7 M Q for a dark matter 
halo consistent with CDM predictions. For comparison, our 
best-fitting NFW model gives M(<0.6 kpc)= 2.0 ±0.6 x 10 7 
Mq, consistent within 2a. Our best-fitting cored profiles in- 
stead give M(<0.6 kpc)= 0.9 ± 0.2 x 10 7 M (constant /3) 
and M(<0.6 kp c) = 1.0 ±0.2 x 10 7 M B , both very similar to 
the estimate of IStrigari et al.l (|2007l ). Therefore in this case 
Sextans would fit in the picture of a minimum mass scale 
for dSphs. Note that the mass within such a small distance 
from the centre, while it appears to come out naturally from 
CDM models, is not very informative on the total dark mat- 
ter content of the satellite: while the enclosed mass within 
0.6 kpc spans 1 order of magnitude (~ 10 7 and ~ 10 s Mq), 
the total dark matter mass may span about 3 orders of mag- 
nitude, as a consequence of the lack of a tight correlation 
between concentration and virial mass, and of the particu- 
lar (a ccretion and orbital) history of the satellite (|Li et al.l 
120091). 

IWalker etahl (2009, 2010) have expanded the above re- 
sults to the mass enclosed within the half-light radius and 
also propose that all dSphs may be embedded in a dark 
matter halo of similar mass (and mass profil e). The mass 
within the half-light radius (estimate from IWalker et al.l 
|2010| . adjusted for the distance to Sextans here assumed) 
that we derive for the best-fitting NFW model (c= 10) is 
M(<0.7 kpc)= 2.6 ± 0.8 x 10 7 Mq, in remarkable agree- 
ment with the determination of I Walker et a l. (2010j3- The 
best-fitting cored dark matter models however give lower 
values, 1.4 ± 0.2 x 10 7 M Q (r c = 3 kpc and /3 = const) and 
1.7±0.4 x 10 7 M Q (r c = 1.5 kpc and /9 = /3 m). These esti- 
mates make Sextan s fall below the "univ ersal" mass relation 
for dSph galaxies in Wal ker et al. I (|2010l ). Until degeneracies 
between different dark matter models are relieved, it is pos- 



8 Also our other NFW models though give a very similar enclosed 
mass (M(<0.7 kpc)= 3.0 X 10 7 M for c= 20 and M(<0.7 kpc)= 
3.2 x 10 7 M Q for c= 35), consequence of the degeneracy between 
concentration and virial mass. 



sible that universal relations between dSphs mass and other 
properties may change according to the adopted models. 

Finally, we would like to point out that for small galac- 
tic systems such as dSphs the total mass is likely to be an 
important factor in determining their evolution. Indeed the 
depth of the potential well will influence how much of the 
gas and metals can be retained and therefore the subsequent 
star formation and chemical enrichement history from the 
system and how much will be ex pelled because of for ex- 
ample supernovae explosions (e.g. Reyaz et al., 20091 ). Since 
outside the last measured point one can only extrapolate 
mass values, it will not be possible to determine total dSph 
masses with the current methods. However, determination 
as far out as possible are still important so as to provide 
constrains for models that aim at reproducing the detailed 
observed properties of stars in dSphs. 



7 SUMMARY 

We obtained VLT/FLAMES intermediate resolution (R ~ 
6500) spectroscopic observations in the NIR CaT region for 
1036 distinct targets along the line-of-sight to Sextans, with 
magnitudes and colours broadly consistent with RGB stars. 
This resulted into 789 stars with S/N and error in veloc- 
ity such as to produce reliable line-of-sight velocities and 
CaT EWs. Among those, 174 are RGB stars probable mem- 
bers of Sextans with line-of-sight velocities accurate to ±2 
km s _1 and CaT [Fe/H] measurements accurate to ± 0.2 
dex. 

The vast majority of the Galactic contaminants could 
be weeded out from the sample using a simple 3-<r kine- 
matic cut. We refined our membership criteria by using the 
Mgl line at 8806.8 A as an empirical indicator of stellar sur- 
face stellar gravity, and therefore distinguish between prob- 
able members of the dSph (RGB stars) and Galactic con- 
taminants (most likely dwarf stars). 

We used the sample of probable members to investigate 
the wide-field metallicity and kinematic properties of Sex- 
tans, and to perform a determination of its mass content. 

Sextans is a metal poor system, with an average 
[Fe/H]= —1.9. Its stars display a wide range of metallicities, 
with the majority being between —3.2 < [Fe/H] < —1.4. A 
spatial variation of the metallicity properties is present, with 
the central parts being in average more metal rich than the 
outer parts. There are indications of a link with the kine- 
matics, with the stars more metal rich than [Fe/H]= —2.2 
having a lower dispersion than the stars with [Fe/H]< —2.4. 

The analysis of the radial trend of the average velocity 
and dispersion of Sextans stars as function of metallicity 
has uncovered the presence of a cold, metal-poor ([Fe/H]~ 
—2.6 from CaT measurements) kinematic substructure at 
projected radii R < 0.22 deg. This structure is consistent 
with being a disrupted stellar clu ster. It is poss i ble th is may 
be the same feature observed bv lKlevna et al.l (|2004l ). 

Our sample of Sextans stars shows a statistically sig- 
nificant velocity gradient along the projected major axis 
of the galaxy and along an axis at P.A.= 191° (8.5 ± 
3.0 km s _1 deg -1 and 7.5^3'q km s _1 deg -1 , respectively). 
Given the current global observational understanding of this 
object it is unclear whether this gradient is due to intrin- 
sic rotation, tidal disruption or geometrical effects. Also, a 



14 G. Battaglia et al. 



better sampling of the outer parts, both in terms of spatial 
coverage and number statistics, would allow to place the am- 
plitude of the velocity gradient on more secure footing and 
perhaps determine the overall velocity pattern throughout 
the face of the galaxy. 

In the hypothesis of dynamical equilibrium, we per- 
formed a mass determination of Sextans, following a Jeans 
analysis and comparing the observed line-of-sight velocity 
dispersion profile to the predictions for cored and cuspy dark 
matter models. Sextans appears to be a very dark mat- 
ter dominated object, with a mass within 2.3 kpc of 2 - 
4x10 s Mg, giving a dynamical mass-to- light ratio between 
460-920 (M/L)v, . Within the explored dark matter mod- 
els and velocity anisotropy profiles, the best fits are given 
by cored profiles with large core radius and mildly radial 
anisotropy. However also a NFW profile with mildly tangen- 
tial anisotropy gives a very good representation of the data; 
the concentration of the favoured NFW is c = 10 somewhat 
lower for what expected for galaxies in the mass range of 
Se xtans when using the virial mass - concentration relation 
of lJind (|2000| ). 
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Parameter 


value 


reference 


I CV TOAAA it 7 f^r\C\r\ P 

y-*j ZvW ) U ,J zUUU / 


IQh 13 m Q 3 s _Q1° 36' 54" 


1 


P.A. 


56° 


2 




0.35 


2 


-encore 


16 6 arcmin 


2 


Rtidal 


160 arcmin 


2 


H. c 


15.5 arcmin 


2 


v sys 


226.0±0.6 km s" 1 


3 


V S ys,GSR 


78.4±0.6 km s^ 1 


3 


(J 


8.4±0.4 km s" 1 


3 




8.4±0.4 km s" 1 


3 


(m-M)o 


19.67 


1 


Distance 


86 kpc 


1 


w 


4.37 ± 1.69 x 10 5 L 


2, 4 


V HB 


20.35 


2 


E(B-V) 


0.0477 


5 



Table 1. The various rows are from top to bottom: coordinates of the optical centre; position angle, defined as the angle between the 
north direction and the major axis of the galaxy measured counter-clockwise; ellipticity, defined as e = 1 — b/a; King core and tidal 
radius; exponential radius; systemic velocity in the heliocentric and in the GSR system; global velocity dispersion, in the heliocentric 
an d in the GSR system; distance modulus and heliocentric distance; lumin osity in V-b and (based on the apparent magnitude measured 
bv llrwin fc Hatzidimitriou 1995 but read justed for a distance of 8 6 kpc bv lLokadE)09T ): V m agnitude of th e hori zontal branch; r edden - 
ing. References: 1 = Mateo (1998); 2 = Irwin & Hatzidimitriou (1995); 3 = this work: 4 = Lokas (2009); 5 = Schlegel ct al. (1998), 
http://irsa.ipac.caltech.edu/cgi-bin/bgServices/nph-bgExec, average value over 5 degrees. 



Field name 


RA(J2000) 


DEC(J2000) 


Date 


Exp. time [sec] 


S/N 


Ntarg 


Sext_100 


10 ft 12 m 58.1 s 


-01°38'05.5" 


2003-12-20 


2700 


18.6 


54 


Sext.101 


10' l l4 m 47.4 s 


-01°11'27.5" 


2003-12-22 


3800 


17.2 


55 


Sext_102 


10 ft 13 m 22.1 s 


-01°21'53.6" 


2004-03-14 


500, 1125, 3600 


2.9,4.5,15.5 


74,92,92 


Sext_103 


10 ,l ll m 13.3 3 


-01°12'01.9" 


2003-12-22 


4000 


13.3 


50 


Sext_104 


10 /l ll m 21.6 s 


-01°31'38.0" 


2004-03-15 


3600 


21.1 


68 


Sext_112 


10 /l 13 m 38.5 s 


-02°00'35.6" 


2008-05-25 


3600 


14.0 


81 


Sext_113 


10 ft 14 m 52.3 s 


-02°24'24.5" 


2004-01-01 


3600 


25.1 


48 


Sext_114 


10 ft 12 m 24.1 s 


-01°58'43.5" 


2004-03-16 


4600 


17.5 


83 


Sext_115 


10 ,l ll m 30.4 s 


-02°09'45.0" 


2003-12-20 


3600 


8.8 


57 


sext_ext_02 


10 /l 17 m 03.1 s 


-01°01'49.8" 


2004-03-17 


3600 


26.8 


76 


sext_ext_03 


10 h 18 m 36.3 a 


-00°24'21.9" 


2004-03-18 


3600 


28.7 


69 


sext_ext_04 


10 ,l 20 m 43.2 3 


-00°06'51.4" 


2004-03-19 


3600 


37.4 


79 


sext_ext_06 


10 ft 13 m 22.9 s 


-00°26'51.1" 


2004-03-20 


4200 


42.0 


68 


sext_ext_08 


10 h 08 m 27.4 s 


-00°52'01.6" 


2004-03-17 


3600 


31.6 


58 


sext_ext_14 


10^07™ 42.2 s 


-03°13'27.8" 


2004-03-20/21 


2511, 3600 


29.9,33.3 


66,66 


sext_ext_16 


10 h 07 m 29 2 s 


-01°40'19.7" 


2004-03-21 


4500 


27.7 


57 



Table 2. Journal of our VLT/FLAMES observations at R~6500 of Sextans. The columns give respectively: the name of the observed 
field, the coordinates of the FLAMES pointing, the UT date of observation, the exposure time, the median signal-to-noise and the number 
of the observed spectra. All these observations were taken in service mode by ESO staff. Field 102 and ext_14 have repeated exposures 
because the first exposures did not meet the required criteria for the maximum allowed seeing and/or airmass. This repeated exposures 
also allow a check for the internal velocity and EWs errors. 







x 2 




P/3 


M(< fl last ) 




M v 


Cored r c 


= 3.0 kpc /3 = const 


0.73 


0.2 


0.06 < p < 0.6 


4.0 ± 0.7 x 10 s 


M 




Cored r c 


= 1.5 kpc p = Pom 


2.5 


0.6 


r a > 1.5 kpc 


3.2 ± 0.7 x 10 s 


M 




Cuspy C : 


= 10 ft = const 


2.9 


0.7 


-1.4 < < 


1.9 ±0.6 x 10 8 


M 


2.6 ± 0.8 x 10 9 M 



Table 3. Parameters of the best-fitting dark matter models for mass modelling of Sextans, when using the observed l.o.s. velocity 
dispersion profile from "Binning2". The columns show: the x 2 i the reduced \ 2 (with the number of degrees of freedom being u = 4), 
the parameter defining the anisotropy (i.e. P itself for the ft = const case, and the anisotropy radius r a [kpc] for the P = /3om case), the 
mass contained within the last measured point (at ~2.3 kpc, assuming a distance to Sextans of 86 kpc), and for the cuspy profile the 
virial mass. 
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Name a. 72 ooo <5j2000 S/N v^Ikms- 1 ] YW [A] [Fe/H] 



sext.000 


10 12 40.27 


-1 29 7.7 


13.0 


241.51 ± 4.64 


2.28 ± 0.46 


-2 34+ ' 27 
z -° -0.31 


sext.OOl 


10 12 47.98 


-1 29 38.8 


22.7 


222.51 ± 1.70 


3.44 ± 0.26 


-1 7S +0.13 


sext.002 


10 12 39.03 


-1 29 59.0 


34.6 


235.84 ± 1.22 


3.81 ± 0.17 


-1 s7 +0.08 


sext_003 


10 13 9.94 


-1 29 34.9 


21.7 


241.29 ± 2.81 


2.93 ± 0.28 


-1 95+ 015 
± - ao -0.15 


sext_004 


10 13 18.97 


-1 26 58.2 


20.0 


216.48 ± 2.40 


1.96 ± 0.30 


9 t-o+0.20 
z.jo_ 23 


sext_005 


10 13 17.14 


-1 26 38.1 


22.9 


245.40 ± 1.48 


2.54 ± 0.26 


9 97 +0.14 



Table 4. This table lists the relevant data for the stars in Sextans observed with VLT/FLAMES which passed our membership criteria. 
The columns indicates: (1) the star ID; (2), (3) star coordinates (right ascension in hours and declination in degrees); (4) S/N/A ; (5) 
heliocentric velocity and its e rror; (6) summed CaT eq uivalent width (EW2 + EW3) and its error; (7) [Fe/H] value and error, derived 
adopting the calibration from [starkcnbur g et al. I <201(t) . 




-40 -20 20 40 -10 -5 5 10 

AV [km s" 1 ] AV/V(a, 2 + ct 2 2 ) 

Figure 1. Comparison between velocity measurements for stars with double measurements in the Sextans dSph. Panel a) Distribution 
of velocity differences for all the stars (dashed line), and for the stars with S/N per A > 10 and estimated error in velocity < 5 km s _1 for 
each measurement (solid line) . The weighted mean velocity, rms dispersion and scaled median absolute deviation fro m the median (where 
the scaled median absolute deviation from the median (MAD) is 1.48xMAD = a robust rms e.g. Hoaglin et al. Il983h are: 0.9 i 3.0 
km s" 1 , 29.0 ± 2.2 km s" 1 , 4.3 km s" 1 (dashed line) and -0.6 ± 0.5 km s^ 1 , 2.5 ± 0.5 km s" 1 , 2.0 km s^ 1 (solid line). Panel b) As 
above but now the velocity difference is normalised by the predicted error. With these S/N and velocity error cuts the measured error 
in the velocity distribution is close to the expected unit variance Gaussian (1.48xMAD = 0.98). 
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Figure 2. Comparison between EW measurements for stars with double measurements in the Sextans dSph. Distribution of summed 
CaT EW (EW2+EW3) differences from Gaussian fit (left) and MglEW differences from integrated flux (right) as a function of S/N 
for the stars with S/N per A > 10 and estimated error in velocity < 5 km s~ -'-for each measurement. Assuming a similar S/N for the 
individual measurements of each star, the error in the summed CaT EW difference is \J2 X oy,w = v2 x 6/ (S/N) while the error in 
the MglEW difference is v2 X oy.w = v2 x 2.8/ (S/N) . The solid lines indicate the f and 3 a region for the error in the difference 
of summed EW. The weighted mean summed CaT EW difference, rms dispersion in the difference and scaled MAD are 0.03 ± 0.05 A , 
0.22 ± 0.06 A , 0.25 A . The weighted mean in the MglEW difference, rms dispersion and scaled median absolute deviation from the 
median arc 0.01 ± 0.02 A , 0.09 ± 0.03 A , 0.08 A . 
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Figure 3. Distribution of heliocentric velocities for the stars observed with VLT/FLAMES in the direction of Sextans and that met our 
selection criteria. The dotted lines indicate the region used for the 3-cr membership selection. 
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Figure 4. [Fe/H] versus GSR linc-of-sight velocity for the observed VLT/FLAMES targets (crosses and squares show the stars outside 
and within the 3<r membership region, respectively; the vertical dotted lines show the velocity membership region). Note that the 
stars found outside the 3<r membership region, which are probable Galactic contaminants, tend to be preferentially found in the range 
—2.0 < [Fe/H] < — 0.5. This means that Sextans stars at [Fe/H] < — 2.0 are less likely to be contaminated by Galactic stars whose velocity 
falls in the 3<r membership selection region. The diamonds indicate stars with velocity consistent with membership but with Mgl EW too 
large for giant stars (see text). The circle indicates one object with velocity and MglEW consistent with membership but with unusual 
position on the CMD (see text). 
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Figure 5. MglEW vs CaT T,W in bins of V magnitude. The range of magnitudes are indicated on the top left corner of the panels. The 
filled red squares show the stars which are highly probable Sextans members (i.e. with velocities within ±2ct from the systemic velocity 
of Sextans) and the open black squares the stars which are highly likely to be MW contaminants (i.e. with velocities at least 4a away 
from the systemic velocity of Sextans). 
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Figure 6. Colour-magnitude diagram of stellar objects for Sextans as obtained from INT/WFC and ESO/WFI imaging data (dots), 
with overlaid the targets of our VLT/FLAMES spectroscopic observations which passed our quality criteria. The squares indicate the 
probable Sextans members, the diamonds the stars that would have been classified as members from the 3<r clipping but that have a 
too large MglEW (> 0.5 A), the crosses indicate the stars with velocities outside the 3cr range of selection. The star encircled has a 
discrepant location on the CMD with respect to the other Sextans members of all metallicities, and it is the outermost star in our sample 
with a projected radius of more than 2 deg from Sextans centre. We believe that this star is likely to be a contaminant as well. 
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Figure 7. Location of the observed VLT/FLAMES fields (solid circles) and targets in Sextans (squares: probable members; diamonds: 
stars with velocities within the 3<r range of selection but with too large Mg I EW for giant stars; crosses: non-members with velocities 
outside the 3<r range of selectio n). The small blue circle indi cate the star with peculiar location on the CMD. The ellipse shows the 
nominal tidal radius (value from Irwin & Hatzidimitriou 1995). 
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Elliptical Radius [deg] No. 

Figure 8. Left: Metallicity distribution as function of elliptical radius for the probable members of Sextans (squares with error-bars). 
The diamonds with error-bars show those stars whose velocities fall within the 3cr range of membership, but that are likely non-members 
based on their large MglEW (> 0.5 A). The small blue circle indicate the star with peculiar location on the CMD. Since the stars 
represented with the diamonds and the blue circle are likely non-members they will not be considered when deriving properties relative 
to Sextans. The red asterisks connected by a solid line represent a running median over 15 stars (except for the last point, which is over 
9 stars. Right: Metallicity distribution for Sextans members (from the squares in right-handside panel). 
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Figure 9. Metallicity distribution for Sextans probable members at projected radii R<0.8 deg (thick red histogram) and at R>0.8 deg 
(thin black histogram). The solid black thin line shows the Gaussian fit to the distribution at R>0.8 deg; the thick solid red line shows 
the best fit to the distribution at R<0.8 deg for a sum of two Gaussians; the position of the central [Fc/H] and dispersion of one of the 
two Gaussian is kept constant to the values found for the distribution at R>0.8 deg. The dashed lines show the two Gaussians separately. 
The error-bars are Poissonian. 
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Figure 10. Velocity field for probable members to Sextans. The velocities are smoothed with a "median filter", i.e. we associate to each 
star the median velocity of the stars located within a square of side 0.2 deg centred on the considered star; this smoothing is just for 
visual purposes and aims at reducing the velocity variations due to the dispersion of the galaxy. We use the velocities in the GSR frame 
to avoid spurious gradients introduced by the component of the Sun and LSR motions along the l.o.s. to Sextans. The colour bar gives 
the velocity scale. The dot ted ellipses are placed at 0.2, .5, and 0.7 deg to give an idea of the distance scale. The solid ellipse shows the 
nominal tidal radius (from llrwin fe Ha tzidimitriou 1995T). 
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Figure 11. i>gsr f° r probable Sextans members versus distance along axes at P.A.= 56°(major axis), 101°, 146° (minor axis), 191° for 
stars located within ±0.15° from these axes (squares with error bars). The probable non-members within the same spatial region are 
shown as diamonds and crosses with error-bars (stars with velocity consistent with membership but with Mg I EW too large for giant 
stars and stars with velocity outside the 3cr selection region, respectively). The solid line indicates the systemic velocity of Sextans; the 
dotted lines indicate a region in velocity of ±5cr from the systemic; the dashed line shows the best-fitting straight line to the data (see 
the text for the values of the best fit). 
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Figure 12. L.o.s. average velocity (top) and velocity dispersion profile (bottom) in GSR system for probable Sextans members versus 
projected radius. Left: profiles derived keeping the number of stars per bin approximately constant to 35 stars per bin; the profiles 
have been derived both when the velocity gradient is included (asterisks: elliptical binning; diamonds: circular binning) and when it is 
subtracted from the individual velocities (squares: elliptical binning; triangles: circular binning). Right: profiles derived for 3 different 
choices of binning with bins of increasing width with projected radius, in the case when the gradient is subtracted from the individual 
velocities and elliptical binning. The horizontal solid line in the top panel indicates the systemic velocity of Sextans. 
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Figure 13. Kinematic properties for probable Sextans members more metal rich than [Fc/H]= — 2.2 (red filled squares) and more metal 
poor than [Fe/H]= —2.4 (blue open squares). We show the variation of the rotation-subtracted GSR velocity versus elliptical radius (a); 
average rotation-subtracted GSR velocity (b); number of stars (c); velocity dispersion profile using rotation-subtracted GSR velocities 
(d). The solid horizontal line indicates the systemic velocity, the horizontal long-dashed lines show the region used for the 3<r membership 
selection. 
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Figure 14. In this figure we summarize the results of the mass modelling when using the observed l.o.s. velocity dispersion profile derived 
with "Binning2" (left: cored dark matter profile; right: cuspy dark matter profile). From top to bottom we plot the reduced \ 2 values for 
the best fit model as a function of core radius (left) and concentration (right), mass within the last measured point at ~2.3 kpc (left) and 
virial mass (right) both in units of 10 8 Mq , velocity anisotropy in the hypothesis of constant anisotropy with radius, anisotropy radius 
in the hypothesis of velocity anisotropy following an Osipkov-Merritt profile. The black squares and cyan triangles show the results in 
the hypothesis of /3 = const and f} = /3om> respectively. For each of the core radius (concentration) values shown in the top left (right) 
panel, the models are optimized over the anisotropy and mass paramaters. The quoted reduced chi 2 values are therefore the value of the 
X 2 for the best fit model over 4 degrees of freedom. 
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Figure 15. Observed l.o.s. velocity dispersion profile (squares with error-bars) overlaid on the best fit model for a cored profile with /3 
constant (solid black line), a cored profile with /?om (solid cyan line), cuspy profile with /3 constant (dashed line) for "Binning2". 
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APPENDIX: COMPARISON WITH WALKER 
ET AL. (2009) 

We compare our dataset of objects with S/N> 10 and veloc- 
ity err ors < 5 km s -1 to the recent data-set of IWalker et al.l 
(2009) consisting of spectroscopic observations of 947 dis- 
tinct targets along the line-of-sight to Sextans acquired with 
the Michigan/MIKE Fiber System (R=20000-25000). These 
spectra cover the wavelength range 5140-5180 A including 
the magnesium triplet absorption feature which the authors 
use, together with calibrating globular clusters, to provide 
an estimate of [Fe/H] for the individual stars. 

Velocities 

There are 141 overlapping ind ividual stars with ve locity de- 
terminations between our and I Walker et all (|2009h samples, 
including probable members and non-members. In Fig. [16] 
(left) we compare the line-of-sight velocities for the stars 
overlapping between the two samples. 

The agreement between the velocity determinations in 
the two data-sets is good: the weighted average velocity dif- 
ference between the two samples is 1.3±1.3 km s _1 , with a 
scaled MAD of 3.8 km s -1 (corresponding to about an aver- 
age velocity error of 2.5 km s _1 in each of the two datasets) 
and 1.4 for the distribution of velocity differences normal- 
ized by the errors. Overall the agreement is good, except for 
5 stars, which have velocities consistent with membership 
to Sextans in ou r sample but helioc entric velocities below 
20 km s" 1 in the IWalker et all (|2009h sample, which would 
classify them as clearly MW stars: the velocities of these dif- 
fer by about 250-300 km s _1 between the two samples. All 
these stars have very reasonable S/N in our data and visual 
examination of the spectra of these 5 stars with discrepant 
velocities does not reveal any problem with the spectrum 
(i.e. CaT is in general visible enough to allow a reasonable 
velocity estimation with our routine). We find it unlikely 
that these may all be binary stars, but the reason for this 
discrepancy remains unclear. 

Excluding these 5 stars does not significantly alter the 
quality of the comparison between the two samples, bringing 
some improvement in the average velocity difference between 
the two samples (0.41±0.44 km s -1 ), while the scaled MAD 
remains almost unchanged. 

Metallicities 

Since the [Fe/H] calibrations in both studies are based on 
different spectral features, a comparison between the [Fe/H] 
determinations of the individual stars is expected to be 
meaningful only when considering probable Sextans mem- 
bers, and not for the MW contaminants for which the lumi- 
nosity correction for the dwarf s tars is likely t o be in correct 
and different in the two studies. IWalker et all (|2009h derive 
[Fe /H] measurements fr o m a M g index and place them on 
the lCarretta fc Grattonl l| 19971 ) scale (their Eq. 7). The cal- 
ibration of the spectral indexes to a [Fe/H] scale requires a 
correction for the V-magnitude of the stars, therefore the au- 
thors provided [Fe/H] values only for those Sextans stars in 
common betwee n their sample and the photometric study of 
iLee et all (| 2003l). using the V-magnitudes derived in the lat- 
ter work. This would result in 50 probable Sextans members 



The DART survey of the Sextans dSph 

with [Fe/H] estimates in common between the two studies. 
In order to increase the number of stars available for com- 
parison with our sample, we check that the V-magnitudes 
obtained f rom our photome try compares well to the determi- 
nations of I Lee et alj (|2003l ) for the same stars, and then use 
our V-magnitudes for those stars in common between our 
sampl e and Wa l ker's but that did not have determinations 
from iLee et all (|2003f) : this gives an additional 44 probable 
members in common between the two samples. The results 
of the comparison are shown in the right panel of Fig. 1161 
in general the determinations from W09 yield higher val- 
ues of [Fe/H] (more metal rich) and the [Fe/H] distribu- 
tio n appears much narr ower. This effect was also highlighted 
bv lWalker etafl l|2009l ) and the authors themselves advise 
against using the determinations from the Mg index as abso- 
lute values. The authors suggest that the Mg spectral index 
could be used as a relative indicator of metallicity. In Fig. [16] 
(right) we show that the metallicities from the Mg spectral 
index correspond to a narrower range than the more reliable 
CaT measurements, suggesting a limited use of the Mg index 
also for relative measurements. It appears also not possible 
to carry out a reliable calibration between [Fe/H] estimates 
from CaT and from Mg index because of the large scatter 
between the two quantities (Fig. [T7J). 

The fact that the Mg index is not a good indicator of 
[Fe/H] prevents us from combining our sample of metal- 
licities with those of I Walker et alj l|2009l ). However, such a 
merged data-set would be most helpful to considerably in- 
crease the size and coverage of both works. 
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Figure 16. Comparison between overlapping stars between this work and IWalker et al ] l|2009h . The comparison for the velocity deter- 
minations is shown in the left panels (134 stars, both probable members and non-members), and for the [Fe/H] values in the right panels 
(94 stars, only probable members). The top panels show the direct comparison, with the solid line indicating the one-to-one relation; 
the bottom panels show the differences among the determinations as a function of the quantities derived in this work, with the solid line 
indicating a null difference. 
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Figure 17. [Fe/H] values as determined in this work against determinations of the Mg spectral index from IWalker et al ] H5oj ) for 94 
Sextans probable members which overlap between the two samples. Note the large range of [Fe/H] values corresponding to a narrow 
range in spectral index values. 



